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Abstract 

In this paper, we present a method to perform renormalized perturbation calcula- 
tions in gauge theories with chiral fermions. We find it proper to focus directly on the 
Ward-Takahashi identities, relegating dimensional regularization into a supplementary 
and secondary role. We show with the example of the Abelian-Higgs theory how to 
handle amplitudes involving fermions, particularly how to handle the matrix 75. As 
a demonstration of our method of renormalization, we evaluate the radiative correc- 
tions of the triangular anomaly in the Abelian-Higgs theory with chiral fermions and 
with Yukawa couplings. Families of chiral fermions with appropriate quantum numbers 
are introduced so that the sum of their contributions to the lowest-order anomaly is 

equal to zero. The left-handed and right-handed fermion of this theory are assumed 

H + 75) (1 — 75) 

to interact with the gauge field in the mixture (1 + 9) h 9 where 9 is 

a constant. This anomaly amplitude is calculated without any regularization and is 
found to vanish. 
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1. Introduction 

Ever since the early seventies, it has become customary to renormalize quantum gauge 
field theories with dimensional regularization, with the classic method of renormalization on 
the basis of subtractions with the aid of the Ward-Takahashi identities largely ignored. In this 
paper, we shall re-examine some basic issues in the theory of renormalization, particularly 
as applied to quantum gauge field theories with spontaneously broken symmetries and with 
chiral fermions. As we all know, a distinctive feature of a quantum gauge field theory is 
that the number of renormalized parameters often exceeds that of the bare parameters. As a 
general rule, a quantum field theory with an excessive number of renormalized parameters is 
likely to be not renormalizable. Take, for example, the theory of scalar QED. It is well-known 
that, if we follow spinor QED to the letter and introduce only two bare parameters, the bare 
charge and the bare mass of the scalar meson, the corresponding quantum field theory is 
not renormalizable. This is because there are three parameters which must be renormalized, 
the third one being the |0 4 | coupling constant. In order to make scalar electrodynamics 
renormalizable, one additional parameter, the unrenormalized |0 4 | coupling constant, must 
be introduced. It is notable that quantum gauge field theories with a spontaneous broken 
vacuum symmetry are exceptions to the general rule. Take the example of the Abelian- 
Higgs field theory. There are in the theory three bare constants in the boson sector: the 
bare coupling of the gauge meson, and the two bare parameters in the Higgs potential- 
same as in scalar QED. These three bare constants generated more than three renormalized 
parameters: the physical masses of the gauge meson and the Higgs meson, the renormalized 
coupling constants of the gauge meson, and the renormalized 3-point and 4-point coupling 
constants of the scalar meson. Of particular importance is the physical mass of the gauge 
meson. Trivially, the mass of the photon in scalar QED is finite, as it is equal to zero. In 
contrast, the mass of the gauge meson in the Abelion-Higgs field theory is non-zero. In 
order to prove that the Abelian-Higgs theory is renormalizable, one must show that this 
extra renormalized parameter of mass is finite to all orders of perturbation. This can be 
done by showing that the gauge meson is related to the other renormalized parameters. 
Such an explicit non-perturbative relation 1 , as well as others which put contraints on the 
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renormalized parameters, are provided by the Ward-Takahashi identities. 

The contents of the Ward-Takahashi identities in a gauge field theory with spontaneous 
symmetry breaking are complicated and require some care to disentangle. In order to sidestep 
these complications, the method of dimensional regularization has been invented 2,3 . With 
the use of this method, the renormalized perturbation series automatically satisfies the Ward- 
Takahashi identities, and one may be led to believe that there is no more need to pay a great 
deal of attention to the identities which have already been incorporated. 

We consider such a belief misplaced. While the method of dimensional regularization 
is convenient to use in explicit calculations, there are two shortcomings associated with 
it. First, the definition of the matrix 75 for a non-integral dimension is subjective and 
controversial up to now 4 . Therefore, the renormalization theory on the basis of dimensional 
regularization alone remains incomplete. Second, and perhaps more important, the Ward- 
Takahashi identities contain rich and physical implications which should be explored. Thus, 
ignoring the Ward-Takahashi identities under the auspices of dimensional regularization is 
not an unmitigated blessing. 

We believe that the Ward-Takahashi identities are the centerpieces of renormalization. 
Their contents should be extracted and utilized. The method of dimensional regularization, 
on the other hand, should be recognized as it is: a mathematical artifice which is helpful to 
use in some cases — no more and no less. 

2. Subtractions in the Fermion Sector 

While the renormalization procedure for the fermion propagator and the fermion vertex 
functions in QED is well-known, new features appear when chiral fermions enter the picture. 
Unlike QED, a theory for chiral fermions has no parity conservation. By Lorentz covariance, 
this 1PI amplitude is of the form 

5m + ma(p 2 ) + m^ 5 b{p 2 ) + c{p 2 )^ + d{p 2 )p^f 5 (2.1) 

where 

8m = m — Too- 
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The unrenormalized propagator S(p) is therefore given by 

S(p) = i [0(1 - c) - m{\ + a) - (bm + #) 75 ] _1 , (2.2) 

which can be shown to be equal to 

q(n) = ■ ^(1 - c) + m(l + a) - (bm + dj>)^ 

lPJ p 2 (l-c-rf)(l-c + rf)-m 2 (l + a + 6)(l + a-6)- 1 ' 

The invariant amplitudes 6 and d are absent in QED, hence there are, in a quantum 
field theory with chiral fermions, two more invariant functions which must be rendered finite 
by renormalization. That all these functions are finite is due to two additional subtraction 
conditions provided by the hypothesis of the existence of in-fields and out-fields. To see this, 
we note that the matrix elements of S(p) are given by 

Safiip) = J d\x - y)e ip ^ < 0\Ti; a (x)My)\0 > . (2.4) 
By assumption, it has a pole at 



where 



Let 



E = JfP + m 2 



R a p = hm j d 3 (x-y)e-^-® < 0\M*)My)\0 > • (2.5) 

y -»-oo 



As y — > — oo, ij)p(y) turns into the in-field and we get 



where the summation over the two spinor states are not exhibited. Also, af n (p) is the creation 
operator for an incoming dressed fermion, and u(p) is the spinor wavefunction for a fermion 
satisfying 

jni(p) = mu(p). 

Similarly, as x — > oo,ip a (x) turns into the out-field and we get 



< OliM*) - / ^ e ~^/^r < 0\u a (p)a out (p)- 



(2tt)3 V E (P) 



The matrix R therefore satisfies 



v(p)Rv(p) = v(p)Ru(p) = u(p)Rv(p) = 0, (2.6) 
where v{p) is a spinor wavefunction for an antifermion satisfying 

ftv{p) = —mv(p). 

Equation (2.6) is obtained as 

v{p)u{p) = u{p)v{p) = 0. 

The pole of S(p) at p = E comes from the integration over very large positive values of 
(x — y ), as the finite range of (x — y ) cannot contribute a value of infinity. Since 

/ dte l(p °- E)t = — —, (2.7) 
Jo P0-E' K J 

iR is the residue of S(p) at po = E. Requiring that the numerator in (2.3) at po = E satisfies 
(2.6), we get 

c(m 2 ) = -a(m 2 ), (2.8) 
d(m 2 ) = 0, (2.9a) 

and 

6(m 2 ) = 0, (2.96) 

which are the subtraction conditions needed for the divergent amplitudes a, b, c and d. 
From (2.3), we find that, when p 2 « m 2 and p 1 is set to m, 

S(p) « -p^Z^m) (2.10) 
p A — m A 

where 

Z ^ = 1 - c(m 2 ) - 2m 2 [c'(m 2 ) + a'(m 2 )} ' 
Eq.(2.11) is obtained by differentiating the denominator of the right-hand side of (2.3) and 
utilizing (2.8) and (2.9). 
Let us define 

S^(p) = S(p)/Z f (m), (2.12) 



then the renormalized fermion propagator S^ r \p) is approximately equal to 

2mi 



p2 _ m 2 



(2.13) 



when p 2 is approximately equal to m 2 and j> is set to m. 

In a perturbative calculation, the 1PI amplitude for the fermion propagator is linearly 
divergent. As we know, a shift of the momentum variable of a linearly divergent integral 
gives birth to a finite term. Therefore, we may interpret a Feynman integral of the 1PI 
amplitude of the fermion propagator as one of symmetric integration, with an unknown 
additive constant arising from an undetermined amount of shift of the momentum variables. 
The symmetrically integrated amplitude is logarithmically divergent. 

There are two observations: (i) Because the amplitudes a and b in (2.1) are multiplied 
by a factor of m, these amplitudes are not linearly divergent by power counting. Therefore, 
they do not contain unknown additive constants, which appear only in the amplitudes c 
and d, (ii) There are no counter terms in the Lagrangian for the amplitudes d and b. Since 
the amplitude b does not contain an unknown additive constant, it must be finite and must 
satisfy (2.9b) on its own. The amplitude d, on the other hand, is allowed an additive constant 
contributed by linearly divergent integrals. This constant is determined by (2.9a). 

The renormalized perturbation series for the renormalized fermion self-energy 1PI am- 
plitude can therefore be obtained as follows. For a graph which has no divergent subgraphs, 
we employ the Feynman rules to obtain perturbatively the amplitudes a, b, c, and d with 
the coupling constants and masses being the renomalized ones. The divergent integrals 
are symmetrically integrated after Feynman parameters are introduced. The amplitude 
ma{p 2 ) + c(p 2 )]/) remains to be logarithmically divergent, and is replaced by the subtracted 
amplitude 

m[a(p 2 ) - a{m 2 )} + [c(p 2 ) - c(m 2 )}^> - 2m 2 [a(m 2 ) + c{m 2 )]{j) - m). (2.14) 

The last term of (2.14) is obtained by requiring that the residue of at p 2 = m 2 with 
p 1 = m is 2mi. These subtractions are the same as the ones in QED. The amplitude b(p 2 ) 
requires no subtraction, while the amplitude d{p 2 ) is replaced by the subtracted amplitude 

d(p 2 )-d(m 2 ), 



so that (2.9a) is satisfied. For graphs with divergent subgraphs, we use the BPHZ formalism 5 . 
We shall express the renomalized 1PI amplitude for the fermion propagator as 

ma r (p 2 ) + m>y 5 b r (p 2 ) + c r (p 2 )^> + d r (p 2 )fty 5 . (2.15a) 

These renormalized functions are related to the unrenormalized functions by 

1 + a r = Z^(m)(l + a), 
l-c r = Z^(m)(l- c), 
b r = Z^(m)b, 
d r = Z^(m)d. 

We shall demonstrate our method of renormalization by the specific example of the 
Abelian Higgs model. The Langrangian density in the Abelian-Higgs model is given by 



(2.156) 



£ = -^F^ + (2V) + W) + E 

i 



+ ^0+0-A o (0 + 0) 2 , 



(2.16) 



with 



= (dp + igoV^tj), 

i>Li = ^(1 + 75)^, 

^Ri = ^(l-Ts)^, 

and 

Ffo, = d,V u - d v Vp. 

In the above, V^, 0, and ipi are the gauge field, the complex scalar field, and the i th chiral 
fermion field, respectively, and 9, g , f 0i , /i and A are constants. The subscript of the 
constants in (2.16) signifies that these constants are bare. The summation in (2.16) is over 
all of the chiral fermions we have introduced into the theory. For simplicity, we shall drop 
the subscript % for the individual chiral fermions from here onward but keep in mind that 
/o, ip, etc refer to any one of the chiral fermions introduced. As usual, we shall put 

v + H + i<p2 
*= V2 
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where v = /xo/\Ao- We shall also call the bare mass of the gauge meson as M = goV 0: and 
the bare mass of the fermion as m = f vo. 

We shall add to the Lagrangian (2.16) a gauge-fixing term and a ghost term. Thus, we 
consider the effective Lagrangian 

L eff = C-^ 2 - iid^Wt) + »M 2 ^ + iag M r){;H, (2.17) 

where 

i = djf* - aM (j)2, 

a is a constant, and £ and r\ are ghost fields. The effective Lagrangian is invariant under the 
following BRST variations: 

8V„ = d^, 5H = g £<f) 2 , 8(f) 2 = -g £(vo + H), 

= 5£ = 0, 
a 

= -igo(l + 0)& L , S^ r = -igoO^R . (2.18) 

Next we turn to the Ward-Takahashi identity for three-point functions involving fermions. 
Starting from 1 

< 0\5Tir](x)4>(y)$(z)\0 >= 0, 
where ip(x) is any of the chiral fermion fields. We get 



'-d^ - M <h)W WIO >-< 0\(irj)(-ig o ^ + 29 + 15 
a 2 



-< 0\(iv)(m9o^ 1 + 26 2 75 )10 > = 0. (2.19) 

We shall adopt the Landau gauge, taking the limit a — > 0. Multiplying (2.15) by e ~ tAx + l P 'v- % v z 
and integrating, we get: 

-iA,r^(p',p) - m ZT^(p',p) = S-\p') l + 29 2 + 1 " - 1 + ^~ 75 ^- 1 ( P ), (2.20) 

where Z = — < 0\(v + H)\0 >. In (2.20), p' and p are the outgoing and incoming 
momentum for the fermion, respectively, and 



The functions T^^, and S(p) are so defined that their lowest-order terms in the 

unrenormalized perturbation series are 7 M — ~—^-,—i / y 5 and - — ■ — , respectively. 

2 p — m 

If we set A = in (2.20), we get 



moZr^wfaP) = -\h^-\p) + S-\p) l5 }. (2.21) 



Let us set, in addition, p 2 = m 2 in (2.21). We get, referring to (2.2), 

m ZT Hlpi p(p,p)\ p 2 =m 2 = -im[l + a(m 2 )]7 5 . (2.22) 

By Lorentz covariance, T^ 2 ^ is a superposition of scalar amplitudes and pseudo-scalar 
amplitudes: 



r ' (j> 2 ii>ip{P ' iP) = 



-t l5 F + F, + t l5 fF 2 + f F 3 + *A(1 + 6)^^G + 



^ 2 ^(m 2 ,m 2 ,0) 
+i£9^G_ + (H + (l + 6>)I±^ + H_9^)(fif - f A) 

(2.23) 

where F iy G±, and H± are invariant amplitudes which depend on p 2 ,p' 2 and A 2 , and P = 
\{p' +p)- Also, F (m 2 ,m 2 ,0) = 1 by definition. By power counting, F 1 is logarithmically 
divergent, while F 2 ,F 3 ,G± and H± are ultraviolet finite. Substituting (2.23) into (2.22), we 
get 

m n 7, 

m[l + a(m 2 )], (2.24) 



Z^(m 2 ,m 2 ,0) 
and 

Fi(m 2 ,m 2 ,0) = 0, z = 1,2,3. (2.25) 

Equation (2.25) with i = 1 insures that the amplitude Fi is actually ultraviolet conver- 
gent. Let us define the renormalized coupling constant / by 



° '^ 2 ^(m 2 ,m 2 ,0)' 

with —= the bare vacuum expectation value of 0. Then we obtain from (2.24) and (2.26) 
v2 

that 

771 

f=-[l + a(m 2 )]Z^m 2 ), (2.27) 
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where the renormalized vacuum expectation value v is given by 



i 



V = V Z/ y/Zfafi). 



By (2.15b), eq. (2.27) can be written as 



i 



/ = -[l + a r (m 2 )]. 



(2.28) 



Thus the values of m and / are related. Hence the renormalized Yukawa coupling constant 
/ is finite if m is finite, and vice versa. 

Let us next study the Ward-Takahashi identity (2.20) in the limit where A is infinitesimal 
but not zero. We shall keep track of terms in this identity which are linear in A. We remark 



as <j>2 is massless in the Landau gauge. A more complete discussion on this can be found in 
Appendix B. We shall only mention here that such a divergence is superficial, as all infrared 
divergent terms in a physical scattering amplitude must cancel. The point is that there is 
no infrared divergence for the physical amplitudes evaluated in the a-gauge with a / 0. 
Since the on-shell amplitudes are a independent 6 , the on-shell amplitudes in the Landau 
gauge are infrared finite. More precisely, at p 2 = p' 2 = m 2 , u(p')T v ^u(p) is infrared finite. 
Nevertheless, some of the matrix elements of T v ^ do have infrared divergence when p 2 and 
p' 2 are both equal to m 2 . Thus the point of subtraction for T v ^ will be chosen to be at 
p 2 = p' 2 = Q 2 , where fl 2 is not equal to m 2 . The physical amplitudes are fl- independent. By 
Lorentz covariance, , r has terms proportional to 7 M , P M , and A M . We shall ignore terms 
proportional to the latter two in the Ward-Takahashi identity (2.20) where is dotted 

with A.. This is because in the limit in which A is infinitesimal with p 2 = p' 2 = Vt 2 



while A^A M is quadratic in A. Thus we shall replace T v ^ in the Ward-Takahashi identity 



that the point of subtraction for T v ^ requires some care. This is because that in the Landau 
gauge adopted in our formalism, T v ^ is infrared divergent if we set both p 2 and p' 2 to m 2 , 




by 



1 



(a+ + 2p + rmi + 0) 



1+75 

2 



1 - 75" 



2 



(2.29) 
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By power counting, the invariant amplitudes (3 + and /3_ are ultraviolet finite. We shall 

therefore pay attention to Z v ^ and a± only. In the lowest order, d + and a„ are both 

equal to 1. We shall define a + to be unity at the subtraction point. Obviously, <5_ is not 

necessarily equal to unity at the point of subtraction. Substituting (2.29) and (2.23) into 

1 + 75 

(2.20) and equating the terms proportional to — + 6) — - — , we get 

Zr{ " !) - l-c r (Q 2 )-d r (Q 2 )-m[l + a r (m 2 )]G + (n 2 ,n 2 ,0). (2.30) 



Eq. (2.30) shows that Z^/Zy^ is a finite number - the counterpart of Z 2 /Zi = 1 in QED. By 
equating coefficients of — — - — in (2.20), in the limit A infinitesimal and p 2 = p' 2 = Q 2 , 

we get 



— -«_(fi 2 , n 2 , o) = i - c r (n 2 ) + d r (n 2 ) - m [i + a r (m 2 )]G4n 2 , n 2 , o), (2.31) 

which shows that a_(fi 2 ,fi 2 ,0) is finite. This also means that the bare parameter 6 needs 
no renormalization. 

Using the rules of subtractions we have given in this section together with the BPHZ 
formalism, one may perform renormalized perturbative calculations to all orders. The renor- 
malized parameters in the Feynman rules are required to obey relations such as (2.30). The 
presence of 75, the meaning of which is controversial in dimensional regularization, presents 
no difficulty whatsoever in our approach. 

One may also make use of the Ward-Takahashi identities for Green's functions of bosons 
to show that the physical constants associated with bosons are finite 1 . Generalization to 
other gauge field theories such as that of SU(2) x U(l) x SU(3) is straightforward. As 
a demonstration, we shall evaluate the lowest-order as well as the next order triangular 
anomalies in the Abelian Higgs theory in the sections below. 

3. The Triangular Anomaly 

A. General Considerations 

In this section we consider the amplitude with three external gauge mesons with polar- 
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izations //, u, p and outgoing momenta k±, &2, ks, respectively, where 

h + h + k 3 = o. (3.1) 

We first observe that this amplitude is linearly-divergent by power counting. We shall 
show, however, that it is not only ultraviolet finite but also has a definite value as a con- 
sequence of Bose statistics. Indeed, there is no undetermined constant in this amplitude, 
despite the appearance of linear divergence. This is true to all perturbative orders. 

To see this, take, for example, the lowest-order term of this amplitude, the diagrams for 
which are illustrated in Fig. 1. We make two remarks: 

(i) Only two of the six diagrams shown in Fig. 1 are topologically independent. However, 
since the amplitudes are linearly divergent, two diagrams with different designation of loop 
momenta yield different amplitudes. Let us denote as M(l, 2, 3) the amplitude corresponding 
to diagram 1(a), where 1, for example, denotes the first gauge meson, with momentum k\ 
and polarization fj,. Then the amplitude in diagram 1(e), with the explicit designation of 
the loop momentum given in the figure, is M(2, 1,3), i.e., it is obtained from M(l,2,3) 
by interchanging (ki, /i) with (A^,^)- The amplitudes for the other diagrams in Fig. 1 are 
obtained similarly. The lowest-order V—V—V amplitude is therefore equal to 

i[M(l, 2, 3) + M(2, 3, 1) + M(3, 1, 2) + M(l, 3, 2) + M(2, 1, 3) + M(3, 2, 1)], (3.2) 

The factor 1/3 is combinatoric, accounting for the overcounting of independent diagrams. 
We emphasize that the amplitude in (3.2) satisfies Bose statistics explicitly, i.e., while k±, k2 
and k 3 are linearly dependent, Bose symmetry is obeyed without the need to invoke (3.1). 

(ii) While the amplitude M(l, 2, 3) is linearly divergent, it is straightforward to show 
that the sum of M(l, 2, 3) +M(2, 1, 3) is only logarithmically divergent. Thus the expression 
in (3.2) is logarithmically divergent by power counting. 

(iii) In order to incorporate Furry 's theorem explicitly, we add to diagram 1(a) the dia- 
gram obtained from it by reversing the direction of the loop. The amplitude M(l, 2, 3) will 
be taken to be half the sum of the amplitudes corresponding to these two diagrams. 
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With these considerations we shall put 

Kvv = e^ihrW^kj, kl kl) + k 2a W\{k\, k\, k\) + k^W^kl kf, k 2 2 )] 

+e^'k la k 2a ,[KW 2 (kl k 2 2 , k 2 3 ) + k p 2 W 3 (kl kl kl) + k p 3 W 4 (kl k 2 2 , k 2 3 )] 
+e^'k 2a k 3a ,[k^W 2 (kl kl k\) + k%W 3 (kl k 2 3 , kf) + k?W 4 (kl kl kl)] 
+^'h a k la ,[klW 2 {kl kl kl) + k^W 3 (kl kl k 2 2 ) + k»W A {kl kl kl)]. 

(3.3) 

This form 7 is obtained by taking advantage of the explicit invariance of (3.2) under the 
cyclic permutation of the external gauge mesons, i.e., 1— >-2, 2— >-3,3— >1. We note that 
all terms in the amplitude are accompanied by the Levi-Civitd symbol e^ upa , as the other 
invariant amplitudes vanish by an extension of Furry's theorem. 

Complete symmetry under all interchanges of the external gauge mesons requires, in 
addition, to require that Tyy V is invariant under 1 44> 2, which leads to 

W 1 (k 2 1 , kl kj) = -W^kl kl kl), (3.4) 

W 2 (kl,klk 2 3 ) = W 3 (kl,klk 2 3 ), (3.5) 

and 

W 4 (kl kl k 2 3 ) = W 4 (kl kl k 2 3 ). (3.6) 

The explicit invariance of (3.3) under the other exchanges of external gauge mesons yields 
no new relations. 

We shall pay particular attention to (3.4), which gives 

Wi(A£0,0) = 0. (3.7) 

In particular, (3.7) gives 

Wi(0,0,0) = 0. (3.8) 

Equation (3.8) together with (3.3) imply that ^yyv together with its partial derivatives 
vanish at k\ = k 2 = k 3 = 0. Since the amplitude Tyy V is linearly divergent by power 
counting, the two conditions of subtraction at k x = k 2 = k 3 = make it finite. As a 
consequence, the evaluation of the VW amplitude requires no regularization. In practice, we 
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integrate over p symmetrically in exactly the same way for all of the amplitudes corresponding 
to the diagrams in Fig. 1. 

The relation (3.8) holds not only for the lowest-order, but in arbitrary orders. This means 
that the amplitude Tyy V is ultraviolet finite to all orders. We note that the Lagrangian does 
not provide a counter term of the form of V 3 . But the amplitude ^vvv needs no counter 
term. Therefore, the Abelian-Higgs theory with a chiral-fermion field is renomalizable. The 
existence of anomalies, which we shall demonstrate in the rest of this paper, is of no relevance 
to renormalizability. 

B. Triangular Anomaly in the Lowest Order 

Over a quarter of a century ago, Adler 8 , Bell, and Jackiw 9 discovered the axial anomaly 10 . 
They considered the amplitude for the triangular graph with two vector vertices and one 
axial- vector vertex. While this amplitude is linearly divergent, enforcement of current con- 
servation on the vector vertices fixes the arbitrary constant of the linearly-divergent integrals. 
They then showed that the current conservation on the axial- vector vertex, argued to be valid 
in the limit m — > as a consequence of chiral invariance, is not satisfied by this amplitude. 

Let us examine the considerations above in the context of the Abelian-Higgs model, in 
which the vector meson couples to the chiral fermion with the form given in (2.16). In this 
model, there is no conservation law either for the vector current or the axial vector current. 
Indeed, since the gauge vector meson couples to the chiral fermion with the form given in 
(2.16), neither the vector current nor the axial- vector current is meaningful individually. 
Furthermore, the quantization of a gauge field requires adding gauge-fixing terms and ghost 
terms to the Lagrangian. The equations of motion containing the contribution of such 
terms are different from the corresponding ones used in the classical version of the theory. 
The current-conservation laws are replaced by the Ward-Takahashi identities, which are the 
consequences of the invariance of the Lagrangian under BRS transformations. The Ward- 
Takahashi identity for the V—V—V amplitude in the Landau gauge is 

ik llM T^ v (ki, k 2 , h) - ZmoF^^h, k 2 , fc 3 ) = 0, (3.9) 

where TyvviX^vv) is defined with a factor <7o(/o<7o) taken away from the V—V—V(<f>T-V—V) 
amplitude. The divergent amplitude ^vvv * s uniquely determined by the requirement of 
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Bose statistics instead of that of vector current conservation. As it turns out, the Ward- 
Takahashi identity (3.9) is not necessarily valid. In the Abelian Higgs model with one chiral 
fermion, an additional term must be amended to the right-hand side of (3.9). This term is 
the ABJ anomaly in the Abelian Higgs model. 

We first calculate the lowest-order term for T^ vv contributed by a single chiral fermion, 
the two Feynman diagrams for which are illustrated in Fig. 2. The lowest-order amplitude 



for TZvvv is e( l ual to 



T 1/p (k 1 , k 2 , k 3 ) = J 



d 4 p 



N a + N b 



(2tt) 4 {p 2 - m 2 )[{p - k 2 ) 2 - m 2 ][{p + k Y ) 2 - m 2 } ' 



where 

N a = Tr 

and 

N b = Tr 



75 (|> + f i + m )r o VP ~h + m )i o VP + m ) 



■ . . v 1 + 29 + 75 . . ,, , . 1 + 29 + 75 . 
75 (-^ + 771)7" -{-j) + h + m )l P s -(-^ - #i + m ) 



2 v r ' ' ' 2 
After some algebra, we get the lowest-order term for T J^yy as 

m r/ . „ Ql ,,„„„/, , r d 3 a5(l — Ea)(l — a± 



(3.10) 



(3.11) 



(3.12) 



:i+9f-9 3 ]e^'k la k 2a , [ 

J m z — 



(3.13) 



where we have symmetrized with respect to ct2, «3 and &2, £3. 

Next we turn to the lowest-order term in the renormalized perturbation series for Tyy V . 
It is convenient to define 



We then have 



L = i(l + 75 ) and i2=I(l- 75 ). 



L 2 = L, 

}2 



R — R, 
LR = RL = 0, 



and 



[l + 29 + l5 ) 



[1 + 6)L + 6R. 



(3.14) 

(3.15) 
(3.16) 
(3.17) 

(3.18) 
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Note that L is the projection operator for a left-handed fermion and similarly for R. A left- 
handed fermion remains to be left-handed after interacting with a gauge vector meson, and 
so does the right-handed fermion. Therefore, a fermion in the loop of the diagrams in Fig. 1 
is either left-handed or right-handed throughout the loop. Consequently, we may calculate 
the contributions from the left-handed fermion and the right-handed fermion separately. The 
contributions from the left-handed fermion is equal to that from a pure left-handed fermion 
(9 = 0) multiplied by a factor of (1 + 9) 3 . This is because there are three gauge vector meson 
vertices on the fermion loop, with each vertex associated with a factor of (1 + 9). Similarly, 
the contributions from the right-handed fermion is equal to that from a pure right-handed 
fermion multiplied by a factor of 9 3 . Consequently, the term in T^ v corresponding to 
diagram (la) is given by 

-Im+flV-flSi f d ' p Tr[ l5l ^ + h)Y(^-hh^] mq , 

2 K + ' l J (2n) 4 [(p-k 2 ) 2 -m 2 }[(p + k 1 ) 2 -m 2 ](p 2 -m 2 y { ' 

In obtaining (3.19), we have averaged over the two directions of the loop momentum. It is 

only after this is done that the argument of the trace in (3.19) has exactly one factor of 75. 

The evaluation of (3.19) is straightforward. We find that the lowest-order amplitude is 

in the form (3.3), with 

W f\klklkl) = -^ [—. [( \t g } 3 ~y d^l-y-a) [< = 1,2,3,4] 

1 v 2 ' 3; 12vr 2 7 m 2 - a 2 a 3 k 2 - a 3 a^ - a^k 2 v ^ 1 L 

(3.20) 

with 

Ni = 2(a 2 - a 3 )a 2 a 3 k 2 + a 3 (3aia 2 + «i«3 + «| + a 2&?)kl 

— a 2 (3a 3 ai + aia 2 + a 2 a 3 + af)kl, 
N 2 = —2a 2 a 3 , 
N 3 = -2aia 3 , 
N4 = 4aia 2 , 

where wf^ is the lowest-order amplitude for Wi defined in (3.3). From (3.20) and (3.13), it 
is also easy to find that the right side of (3.9) should be amended with an anomalous term, 
the lowest-order of which is 

e ^-\^ 2CT ,_L[(i + of - e% (3.21) 

16 



in the Abelian Higgs theory with one chiral fermion field. 



4. Anomaly of the Next Order 

A. General Considerations 

It is customary to introduce a family of chiral fermions with appropriate quantum num- 
bers in such a way that the sum of their contributions to the lowest-order anomaly vanishes. 
With this being done, let us go on to study the anomaly in the next order. 

In the next two sections of this paper, we shall calculate the radiative corrections of the 
triangular anomaly in the Abelian Higgs theory from a single chiral fermion. The contribu- 
tions from the family of chiral fermions can be obtained by adding the contributions from 
each of the chiral fermions. As we shall see, this former anomaly, expressed in terms of 
renormalized parameters, is equal to zero 12 . 

Before we plunge into the details of the calculations, a few words of perspective may be 
in order. The anomaly amplitude is a Feynman integral. Because of the symmetry of the La- 
grangian, the integrand for the anomaly amplitude consists of terms which, were we allowed 
to change freely the designation of loop momenta for each of them, would cancel completely. 
However, the anomaly amplitude is a linearly divergent integral, and shifting the variable of 
integration of such an integral without compensation is not legitimate. We shall call such 
illegitimate cancellations superficial cancellations. In the lowest-order anomaly amplitude, 
there is only one loop-momentum, and the amount of shifts required for cancellation is a 
linear superposition of the external momenta. Consequently, the invariant amplitude of the 
lowest-order anomaly can be expressed as an integral over a surface on which the loop mo- 
mentum is infinite, as commonly done in the literature. As a result, the anomaly amplitude 
is a constant, as the external momenta which appear in the integrand of this amplitude can 
be dropped. In the anomaly amplitude of the next lowest-order, there are two loop-momenta, 
and the shift for one of the loop-momentum required for cancellation may involve the other 
loop-momentum. Nevertheless, since the Jacobian of integration for such a shift is equal to 
unity, the contributions to the anomaly amplitude again come from a region in which the 
loop momenta are infinitely large. In this region, the external momenta are again negligible, 
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hence the anomaly amplitude in the next lowest-order is also a constant. 

B. Reduction of Amplitude 

The relevant diagrams for Tyy V in the lowest and next to the lowest orders are illustrated 
in Fig. 3 and Fig. 4. The diagrams in Fig. 3 are just the lowest-order triangular diagrams 
with insertions of radiative corrections to the vertices and internal lines. Instead of using six 
diagrams of this kind with permuted external lines, (and dividing the sum of the correspond- 
ing amplitudes by 3), we shall only use two diagrams here. Bose symmetry is maintained 
by an ingenious choice of the loop-momenta 13 . Referring to Fig. 3, we see that, with the 
loop-momentum so chosen, diagram (3a) and diagram (3b) are individually invariant under 
the cyclic permutation of 1 — > 2, 2 — > 3, 3 — > 1. In addition, these two diagrams turn into 
each other under 2^3. Therefore, with the understanding that the integration over p is 
symmetrically performed in exactly the same way, the sum of the amplitudes corresponding 
to these two diagrams is completely Bose-symmetric. 

Next we turn to the diagrams in Fig. 4. We note that each of the first two diagrams in 
Fig. 4 has four internal fermion lines, while each of the rest of the diagrams has three internal 
fermion lines. The diagrams to be included in our calculations for Tyy V are obtained from 
diagrams (4a)-(4f) by permuting the external lines in all possible ways. 

The diagrams relevant to the amplitude T^w are obtained from those for the amplitude 
IYvv by replacing an external V line by an external 02 line, plus the diagrams obtained 
from diagram (4g) by permuting the external lines in all possible ways. 

Since there are a large number of diagrams for each of the amplitudes Tyvv and T^w, 
the calculations of these individual amplitudes are tedious. However, it is much simpler to 
calculate the anomaly amplitude contributed by these diagrams, as we may take advantage 
of the Ward-Takahashi identities such as (2.20). 

Let us first consider the anomaly amplitude (the expression we should add to the right- 
hand side of (3.9)) contributed by the diagrams in Fig. 3. More precisely, we calculate T vvv 
contributed by such diagrams and T^yy contributed by those obtained by replacing the 
gauge meson of momentum k\ with a 02 meson in these diagrams. By (2.17), the anomaly 
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amplitude from diagram (3a) is: 
d A p 



+ 



I 



(2n)* 
dp 



Tr 



Tr 



r"(p - k 23 , P - k l2 ) l + 2 i + ^ s{p - k 31 )r p ( P - k 31 , P - k 23 )s( P - k 23 ) 



r"(p - k 23 ,p- k 12 )S{p - k 12 ) 



1 + 20 



T p (p - k 31 ,p - k 23 )S(p - k 



2:5 I 



(4.1) 



and the anomaly amplitude from diagram (3b) is 
-Tr 



d A p 
dp 

In the above 



T U ( P + fcia, p + k 23 )S(p + k 23 )Tf(p + k 23 , p + k 31 ) 1 + 29 n + l5 S(p + k l2 ) 



Tr 



T v {p + k 12 ,p + A; 2 3).S(p + k 23 )T p {p + /c 23 ,p + fc 3 i)-5(p + k 3l ) 



1 + 2^-75 



(4.2) 



fcjj — {ki kj) /3. 

The four terms in (4.1) and (4.2) are schematically represented by the four diagrams in Fig. 
5. We note that in these diagrams, the dotted line located behind (in front of) the vertex 

carries a factor — ^ 5 (— — ), the order being dictated by the direction of the 

arrow. (See the right-hand side of (2.20) for the origin of the factors (l + 26>±7 5 )/2). We also 
mention that the momentum variables associated with a vertex which is joined by a dotted 
line are specified by the momentum of the external gauge meson and the momentum of the 
fermion line on the opposite side of the dotted line. For example, for the vertex function T u 
in diagram (5a), the outgoing momentum is p — k 23 , while the incoming momentum is not 
p — k 3 i, but is 

(p - ^23) + k 2 =p- k 12 . 

Let us change the loop momenta in diagram (5a) and in diagram (5d) so that the momenta 
for both the upper fermion lines are designated as p. Similarly, we change the variable of 
integration for diagrams (5b) and (5c) so that the momenta for both lower fermion lines in 
these diagrams become p. The anomaly amplitude given by the sum of (4.1) and (4.2) then 
becomes 

Ji + J 2 (4.3) 
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where J\ is a volume integral given by 



with 



N R = 



(2k) 



-r"(p + A; 3 ,p - fci)^(l + 20 + 75) + ^(1 + 29 - 75 )r"(p - fc 2 ,p) 



(4.4a) 



5(p)^(p,p + A;3) ( 5(p + A;3), 



(4.46) 



and 

7V 6 = n P ,p + £; 2 )S(p + £; 2 ) 



-r-( P + k 2 ,p - k 1 ) 1 -±^ + I + ^I^ - fc.p) 



5(p). 
(4.4c) 



The term J 2 comes from changing the variable of integration for the linearly divergent inte- 
grals, and is given by the surface integral 

1 r dS" 



J? 



3 J (2tt) 4 



7: 



(4.5a) 



where N 7 is given by 



N 7 = (h - k^Tr^ip + h,p- h) ( 1 + 29 + ^ S(p)r»(p,p + k 3 )S(p + k 3 )\ 



+ (k 2 - k 3 ) fl Tr[ {1 + 26 ) l5) T»(p - k 2 ,p)S(p)TP(p,p + k 3 )S(p + k 3 )} 



(h - k 2 ) lx Tr[T v {p,p + k 2 )S{p + k 2 )r»(p + k 2 ,p- h 



(1 + 20 + 75) 



Sip)] 



(4.56) 



- (k 2 - k 3 ) ll Tr[T»((p,p + k 2 )S(p + k 2 ) {1 + 26 n 7s V (p - k 3 ,p)S(p)], 



and where the integration is over a three-dimensional surface in the p-space at infinity. 

Since we are interested in the anomaly only in the next lowest order, (4.4b) and (4.4c) 
can be replaced by 

(1 + 20 + 75) (1 + 20 -75). 



Vip + k^p-kJ- 



-r»(p-k 2 , P ) 



-r 



,(l + 20 + 75 ) 



(4.6a) 



and 



N 6 = 7 



p 1 — m 
,(1 + 20 + 75) 



]/> + ^ 2 - m 



r^ip + k^p-h) 



(1 + 20 + 75) 



(1 + 20 -75) 



r p (p-k 3 ,p) 



(4.66) 



f — 7X1 
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as the difference of the two vertex functions inside each of the brackets is already at least of 
that order. 

The vertex function T in (4.4), (4.5) and (4.6) is that of the V^^-vertex. We shall express 
this vertex function in a renormalized perturbation series. Since the counter term for this 
vertex function is equal to a constant times 7^- — - ^ , it is easy to prove that the 
counter terms for the T function in (4.6) cancel. 

We draw in Fig. 6 diagrams for the lowest order radiative correction of T v ^. In partic- 
ular, diagram (6a) depicts the radiative correction to T v ^ due to the exchange of a vector 
meson. We draw in Fig. 7 schematic diagrams representing Ji of (4.4), with the diagram (6a) 
inserted for the vertex sub-graph. Note that the dotted lines in Fig. 7 are always joined to 
bare Vipip vertices, each of which is assigned a factor 7 M [(1 + 9)L + 9R] in the Feynman 
rules. It is easy to see that the amplitude corresponding to diagram (7a) is equal to that 
corresponding to diagram (7d) after the change of integration variables 

p — >■ p - q - k 2 

is made. This does not mean, however, that these two divergent amplitudes cancel each 
other, as a change of integration variables is not allowed. Similarly, the anomaly amplitude 
corresponding to diagram (7b) is equal to that corresponding to diagram (7c), after a change 
of variable is made. The sum of the amplitudes corresponding to these four diagrams will 
be evaluated in the next section. 

Next we consider the anomaly amplitudes corresponding to Ji of (4.4) with diagram (6b) 
inserted for the vertex. These anomaly amplitudes are represented by the diagrams in Fig.8. 
We emphasis two points: (1) The integration over q, the loop momentum for the radiative 
correction of the vertex, is to be performed before that over p, the loop momentum in the 
triangular diagrams in Fig. 3. This order of integration must be followed strictly. This is 
because individual Feynman integrals are divergent and the value of an individual integral 
depends on the order of integration. (2) The integration over q for the sum of the anomaly 
amplitudes corresponding to the diagrams in Fig. 9 is convergent. This is because the ultravi- 
olet divergent in the g-integration for diagram (8a) cancels that for diagram (8b). Similarly 
for the ultraviolet divergences in the g-integration for diagrams (8c) and (8d). Therefore, we 
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are allowed to choose to symmetrically integrate over q for all of these amplitudes. 

The anomaly amplitudes represented by the diagrams in Fig. 8 are cancelled by corre- 
sponding terms in the non-planar diagrams of the form of Fig. (4a). To see this, we draw all 
twelve of these non-planar diagrams in Fig.9. Notice that only six of them are topologically 
independent. For example, diagram (9a) is topologically the same as diagram (9h). We draw 
all twelve of them for the sake of symmetry and, to compensate for overcounting, we shall 
multiply the amplitude for each of these diagrams by a factor of 1/2. Notice also that we 
have assigned the loop momenta in these diagrams in such a way that the sum of amplitudes 
from these six diagrams are Bose-symmetric with respect to the exchange of any two of 
external vector mesons. The Feynman integral for each diagram contains two integrations, 
namely q and p. We shall choose to integrate over q first. (The result is independent of the 
order of integration, as long as we follow the same order for all the diagrams in Fig. 9.) We 
also mention that, power counting, the integration over q is convergent for the amplitude 
of an individual diagram. We shall therefore choose to integrate symmetrically over q. The 
subsequent integration over p will also be chosen to be symmetric, as is discussed at the end 
of Sec.3A as well as in Sec. 5. The contribution of diagrams (9a) and (9b) to the anomaly 
amplitude is schematically evaluated in Fig. 10 and Fig. 11. A dot next to the external line 
of vector meson 1 in these diagrams represents the multiplication of iki^ to the VVV ampli- 
tude. The anomaly amplitudes are represented by the diagrams following the equality signs. 
The diagrammatic equality is obtained via the Ward-Takahashi identity 



ai, r Mo) _ 2// 2 r (0) - ? n w ' - i n (U> ' (a 7^ 

where the superscript (0) on a symbol signifies that the quantity the symbol represents is 
that of the lowest order. It is easy to see that the sum of the anomaly amplitudes from 
diagram (10a) and diagram (8b) is equal to the anomaly amplitude from diagram (8b) with 
the 75 associated with the dotted line in diagram (8b) deleted. Similarly, adding the anomaly 
amplitudes from diagrams (10c), (11a) and (11c) simply eliminates the 75 associated with 
the dotted lines of diagrams (8a), (8d) and (8c) respectively. It is also straightforward 
to show that the sum of these anomaly amplitudes is equal to zero. Briefly, the anomaly 
amplitude from diagram (8a), with the aforementioned 75 deleted, is actually a convergent 
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)(0) 



)(0) 



integral. Therefore, it is legitimate to change the designation of the loop momenta and this 
anomaly amplitude cancels the corresponding one from diagram (8d). Similarly, the anomaly 
amplitudes from diagrams (8b) and (8c), with the aforementioned 75 deleted, cancel each 
other. 

We may consider in the same way the anomaly amplitudes from the diagrams obtained 
from those in Fig. 8 with the 2 Tines replaced by the if-lines. These anomaly amplitudes 
are cancelled by the anomaly amplitudes from diagrams (10b), (lOd), (lib) and (lid). 

Next we consider the diagrams in Fig. 5 with the diagrams (6f) and (6g) inserted for the 
vertex. The resulting diagrams are illustrated in Fig. 12. We note that two of the external 
lines attached to the triangular fermion loops in these diagrams are scalar-meson lines. It 
is therefore tempting to argue that, since the amplitude corresponding to such a triangular 
loop vanishes as a consequence of Lorentz invariance, the anomaly amplitude corresponding 
to a diagram with such a loop also vanishes. Such an argument is incorrect, as we do not 
integrate over the momentum of the fermion loop first. 

We draw in Fig. 13 some of the anomaly diagrams derived from the Feynman diagrams 
(9c)-(9f) and (9i)-(91). Note that these anomaly diagrams are topologically equivalent to the 
anomaly diagrams in Fig. 12. The corresponding amplitudes differ, however, as the loop 
momenta in these diagrams differ. We shall show in Appendix A that those diagrams in Fig. 
12 gives zero contribution to the anomaly amplitudes while those in Fig. 13 give nonzero 
contribution. There are also other Feynman diagrams which can be reduced to diagrams as 
those shown in Fig. 12 which we do not show here since they also give zero contribution to 
the anomaly amplitude . The sum of the anomaly amplitudes from the diagrams in Fig. 13 
is equal to 

J 3 = e^'^ff 2 9(l + e). (4.8) 

Next we consider the anomaly amplitude corresponding to the diagrams in Fig. 7 with 
diagram (6d) replacing diagram (6a) as the vertex subgraph. It is straightforward to see 
that these anomaly amplitudes are exactly canceled by the anomaly amplitudes generated 
by diagram (4b) with the vector meson of momentum k\ attached to the fermion loop, and 
diagram (4d) with the vector meson of momentum k\ at VfoH vertex. To demonstrate this, 
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we draw in Fig. 14 a set of such anomaly amplitudes which exactly cancel. Diagram (14a) 
is diagram (7a) with diagram (6d) replacing (6a) as the vertex subgraph. Diagram (14b) 
is an anomaly diagram derived from diagram (4b). Diagram (14c) is an anomaly diagram 
derived from diagram (4d). The factors associated with the dotted lines are - — — — — - ^ , 

— , and (— %){— 175), respectively. The sum of these three terms is zero. 



2 

In summary, the anomaly amplitudes from various diagrams generally cancel as a con- 
sequence of the Ward-Takahashi identities. There are four exceptions. The first exception 
is the contribution from the nonplanar diagrams Figs. (9c), (9f), (9i) and (91) and is given 
by J 3 . The second exception is the volume integral J\ corresponding to the diagrams in 
Fig. 7, for which changing integration variables to effect cancellation is illegitimate, as we 
have discussed. The third exception is the surface integral J2 given by (4.5). J\ and J2 
will be evaluated explicitly in Appendix E. Referring to (4.8), (E.20), (E.30) and (E.31), we 
obtain for the sum of the above three exceptions 14 

■h + h + J 3 = 0. (4.9) 

The last exceptions are the anomaly amplitudes derived from the diagrams (4c) and (4e). In 
each of these diagrams, there is a triangular loop with three vector meson lines attached to 
it. The anomaly amplitudes due to this triangular loop have been calculated in Section 3B. 
This anomaly amplitude vanishes as we sum over the contributions from the family of chiral 
fermions with appropriate quantum numbers. 



5. Summary 

In this paper, we have demonstrated how to obtain renormalized Ward-Takahashi identi- 
ties for the Abelian-Higgs theory with chiral fermions. With the aid of this set of identities, 
one can calculate unambiguously physical quantities by the classic method of renormalization 
on the basis of subtractions. As an example, we calculated the next order triangle anomaly 
in this theory which is shown to vanish 14 . This calculation is performed without the intro- 
duction of regularization. Some of the calculations on higher order triangle anomaly in the 
past are based on the introduction of cutoffs in the Feynman propagators and thus the result 
of such a calculation depends critically on the method of regularization. For example, it is 
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well-known that shifting the loop-momentum of a linearly divergent Feynman integral must 
be compensated by the addition of a finite and nonzero constant. However, if we regularize 
the propagators in this Feynman integral, no such compensation is necessary for such a shift. 
More precisely, 

d 4 p f (p + a)^ A p^ A _l 

(5.1) 



A^Lj (27r) 4 I [(p + a) 2 - m 2 } 2 A - (p + a) 2 (p 2 - m 2 ) 2 A - p 2 
d 4 p f + A Pn A 



/ f dp f 



(2tt) 4 a^oo I [(p + a) 2 - m 2 ] 2 A - (p + a) 2 (p 2 - m 2 ) 2 A - p 2 . 
The left-hand side of the above equation is equal to zero as each of the two integrals on 
the left-hand side is convergent when A is finite. On the other hand, each of the two terms 
in the integrand of the right-hand side, with A set to infinity, is the integrand of a linearly 
divergent integral, and hence the right-hand side is not equal to zero. Generally, caution 
must be exercised when the integrand of an integral depends on a cutoff parameter. The 
limit of such an integral is not always equal to the integral of the limit of the integrand as 
the cutoff parameter is made to go to infinity. This means that the value of a divergent 
integral is dependent on the particular method of regular izat ion. 

In contrast, we start with a rigorous theory on the basis of a gauge invariant Lagrangian 
with the symmetry of the vacuum spontaneously broken. We canonically quantize this 
theory, adding gauge fixing terms and introducing the associated ghost fields. The ghost 
terms are so constructed that the Lagrangian has the BRST invariance. As a consequence 
of the BRST invariance, a system of Ward-Takahashi identities is satisfied by the Green's 
functions of this theory. While the quantum theory so obtained is ultraviolet divergent, the 
Ward-Takahashi identities are taken to form the guiding rules of renormalization. We obtain 
the value of the anomaly amplitude via a rigorous and unambiguous scheme of renormalized 
perturbation. There is no need to make up an axial vector current which is not in the 
theory. In our calculation of the anomaly amplitude, no regularization is needed. While a 
VW-amplitude for an individual Feynman diagram can be linearly divergent, the full Bose 
symmetry with respect to all the three external vector mesons makes the VW-amplitude 
ultraviolet finite without the need of introducing a parameter of regularization. The crucial 
fact is that (3.3) and (3.8) enforce Tyy V as well as its partial derivatives to vanish at the 
point where all external momenta are zero. One can therefore make a subtraction at this 
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point and the anomaly amplitude is then explicitly utlraviolet convergent. 

In practice, we choose the loop momenta for the diagrams in such a way that full Bosc 
symmetry is explicitly satisfied, as was illustrated, for example, in the diagrams in Fig. 3. 
We obtain the anomaly amplitude by integrating over such loop momenta with subtractions. 
Since the subtracted integral is ultraviolet finite, it is mathematically rigorous to adopt 
symmetry integration for the subtracted integral under which both the subtracted terms 
and the unsubtracted terms turn out to be finite. Because the subtracted term satisfies 
the full Bose statistics explicitly on its own, it is equal to zero. The anomaly amplitude is 
therefore equal to the unsubtracted term integrated symmetrically with the loop momenta 
chosen as in Fig. 3 and in this case is found to vanish. 

The renormalization procedure in the above can easily be extended to nonabelian Higgs 
theories with chiral fermions such as the standard model. The renormalized Ward-Takahashi 
identities so obtained can be used to obtain exact relations among various physical quantities 
that are of experimental interests. Thus one can check the validity of the standard model 
against precision measurements that have been done or will be done in the near future. 
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Appendix A 



In this Appendix, we calculate the anomaly amplitude corresponding to the diagrams in 
Figs. 12-13. 

All of the anomaly amplitudes from the diagrams in Fig. 12 are equal to zero. To show 
this, let us remind ourselves that we are going to symmetrically integrate over q first and 
then symmetrically integrate over p. To illustrate how this works, it is sufficient to only look 
at diagram (12a). The Feynman integral for (12a) is 

d 4 p d A q M 



where 



i9(i +e)r 



M ={2q + k 2 ) u Tr 



(2tt) 4 (2tt) 4 V ' 



= AieP a ^ Pa k w (q - fc 1 ) ff (2g + k 2 f . 
V = [(q + h) 2 - 2/i V [(p + h) 2 - m 2 } [p 2 - m 2 } [{p + q- k x ) 2 - m 2 



(Al) 



(A2) 



(A3) 



It turns out that the anomaly amplitude given by (A.l) is independent of the masses in 
the theory. To simplify the presentation, we shall drop the mass terms in (A. 3). Next we 
introduce Feynman parameters and reduce (A.l) into 

d A p d 4 q N 



24^(1 + 6)f j d 5 a5{l - ]T a) J 



(2tt) 4 (2tt) 4 D 



where 



D =a 1 (q + k 2 ) 2 + a 2 q 2 + a 3 (p + k 3 ) 2 + a 4 p 2 + a 5 (p + q- k x ) 2 

F(p) 

= («i + a 2 + a 5 ) (q - 5q) 2 + 



where 



5q 



ai + a 2 + a 5 ) ' 
Oib{ki - p) - aik 2 



(ai + a 2 + a 5 ) 
F(p) = r(p - 5p) 2 + (ai + a 2 + a 5 )A, 

5p = -\a 5 (ai + a 2 )ki + 0:10:5^2 - a 3 (ai + a 2 + a 5 )k 3 ], 
r 

r = a 5 (ai + a 2 + a 3 + 0:4) + («i + a 2 )(a 3 + a 4 ), 



1 



A = -\ a 2 a4a5k 1 + \a 2 a 3 a^+aia 2 (a3+a4+a5)]k 2 + [aia4a5+a 3 a4 : (ai+a 2 +a5)]ki 
r I 



(A4) 

(A5) 

(A6) 
(A7) 
(A8) 
(A9) 
(A10) 
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Since the integration over q is convergent, we make a shift in q — * q + 5q and then integrate 
over q. Upon shifting q, the numerator M takes the form 

M = 2ie vap/3 q 2 p a k 3p + Aie pa " a (8q - h^ppk^Sq + k 2 f . {All) 

(A. 4) becomes, after integrating over q 

30(1 + 0)f f , 5 c n _ V x A <* 4 P f ^ | (<*i + a2 + a 5 )N 2 \ 

n2 jaaO{L 2^ a )J (2 7 r)n(ai + a 2 + a5)i ?2 (p) ^ 3 b) " J' 

(A12) 

where 

N 1 = 2i^ Pa k 3p , (A 13) 

AT 2 = 4ieP°P°(6q - k^ppk^d q + fc 2 ) v . (A14) 

Now we write out explicitly 5g in terms of p and make a shift in p — ■> p + 5p. Here, one 
should be careful about shifting p in (A. 13) since the it involves a linearly divergent integral. 
Thus, when we integrate out the first term in (A. 12), it contains two pieces, the first piece is 
a volume integral over p and the second piece is a constant due to shifting p by an amount 
5p. The second piece is equal to 

39(1 + 9)f 2 e uaf,/3 k la k 3/3 f d 5 a 5(1 - £ a) a 2 a 5 



167T 4 J (a\ + a 2 + as) r 

The first piece is equal to (A. 12) with N\ replaced by 



r d°ad{L - £aj ^2^5 (A 15) 



He^^fciaV (A16) 

and -/V 2 replaced by 

-2ie^ -. j a2a " - 2 k la k 3 ^. (A17) 
(ai + a 2 + a 5 ) 2 

Substituting (A. 16) and (A. 17) into (A. 12) will give an integral of the same form but with an 
opposite sign as that of (A. 15) and so the expression in (A. 12) is equal to zero. All the other 
diagrams in Fig. 12 can be handled in a similar fashion and one will get zero contribution 
for each of them. 

Next we look at diagrams (9c)-(9f) and (9i)-(91). It is easy to see that after being 
multiplied on by (9d), (9e), (9j) and (9k) are reduced into diagrams similar to that of 
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Fig. 12. For example, (9d) is reduced to (12a) and (12b), similarly for (9e), (9j) and (9k). 
Therefore, the diagrams (9d), (9e), (9j) and (9k) give zero contribution. The only diagrams 
that would contribute to the anomaly amplitudes are (9c), (9f), (9i) and (91). The anomaly 
amplitude from these diagrams are schematically shown in Fig. 13. For example, (9c) is 
reduced into diagrams (13a) and (13b), etc. Note that even though Fig. 13 has the same 
topology as that of Fig. 12, the momentum assignments are different. For example, anomaly 
amplitude corresponding to diagram (13a) is 



where 



N = (2q + k 2 ) v Tr 



7s7 P (^ + i)1>{jt> + 4- 



(A19) 



D = [(q + k 2 ) 2 - 2/2 2 ]q 2 [(p + q) 2 - m 2 }[(p + q- k 3 ) 2 - m 2 }(p 2 - m 2 ). (A20) 

In the above, we can see that the momentum q appears in four of the five propagators in 
(A. 20) as opposed to three in (A. 3). By restricting ourselves to integrating over q first and 



then over p similar to the above, one gets a nonzero result for this diagrams which is equal 

512tt 4 J v ; 
contribution from Fig. 13 is 



to e vpaa ' f 2 6(l + 6). Since each diagram in Fig. 13 gives identical result, the total 



J 3 = e^'^ff 2 9(l + 9). (A21) 
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Appendix B 

The renormalized vertex function is equal to the unrenormalized vertex function 

multiplied by a renormalization constant Z v ^. Customarily, one chooses the point of sub- 
traction at A = 0, with both p and p' on the mass shell, i.e. p 2 = p' 2 = m 2 . However, in 
the Abelian Higgs theory formulated in the Landau gauge, the 02-fi+ld is massless. As a 
consequence, the vertex function T v ^ is infrared divergent when p and p' are both on the 
mass shell. Take, for example, diagram (6b). The inverse of the two fermion propagators 
corresponding to this diagram have the following factors at p 2 = p 2 = m 2 , (where pf = p' 
and pi = p) 

(Pi + q) 2 -m 2 =2pi-q + q 2 , 

(p f + q) 2 - m 2 = 2p f ■ q + q 2 . 
In addition, the inverse of the 02-propagator is equal to q 2 . As a consequence, the amplitude 
corresponding to diagram (6b) is divergent at q — 0. 

While T v ^ considered as a matrix is infrared divergent, some of its matrix elements are 
infrared finite. In particular, if we insert T v ^ into the two bispinors u(pf) and u(pi), the 
amplitude is infrared finite. As an example, let us consider the amplitude corresponding to 
diagram (6b). The numerator for the amplitude is 

75 0*/ + 4 + m)^ 1 + + 75 (^ + i + m)75- 

Since 

U>i + 4 + m) 75 = 75 ( -4 + m), 
this numerator operating on u(pi) becomes 

ti j \ „ 1 + 20 + 75, ,x 

7 5 (^ + 4 + m)Y — H-i), 

which vanishes at q — 0. Therefore, the matrix element of T v ^ between the bispinors u(pf) 
and u(pi) is infrared finite. As a matter of fact, we also have 

uiPfhs^f + 4 + m) = u(pf)(-4)l/5- 
Therefore, the numerator vanishes as q 2 when q — > 0. 
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Since the vertex function T v ^, considered as a matrix, is infrared divergent, the point of 
subtraction will be chosen to be at A = and p 2 = p' 2 = Q 2 , with Vt ^ m. The renormalized 
vertex function is therefore equal to Z V ^(Q 2 , Q 2 , 0) times the unrenormalized vertex 

function. 



(1) Diagram (6a) 

The diagram (6a) in Fig. 6 gives the following renormalized amplitude to Tf, 



d 4 q 



(2vr) 4 {{p f + q) 2 - m 2 }{{ Pi + q) 2 - m 2 }{q 2 - M 2 ) 
—term of subtraction, {B.l) 



where 

^ 1 + 20 + 75 i + 2 g + 75 a , l + 26 + 7 5 ( q a q al 
n a = 7 (p f + i + m)Y (p~ t + i + 771)7 W - -~r 



By making use of (3.15)-(3.18), we may reduce the numerator n a in (B.2) to 



i • 0r ! i— ^ • (>-> l 



2 2 
+0(1 + 0)m{ 7 V(A + + m)Y' [(1 + 0) 2 



(B.2) 



1 +75 gl -75 



0- 



+ + 



2 

1-75 



(B.3) 



2 ' 2 

We can further reduce (B.3) with the following identities 

Ivl^lula = ^9p,v and 7^7a = -27 M . 

Substituting (B.3) into (B.l) and introducing Feynman parameters, we find that the inte- 
gration over q is convergent. We get 



->/«(<*)_ 



a ao(l — a) h / a ad (I — a) I — — H — - 

c a J y a a a a 



0(1 + 0)m{ / rf 3 a5(l - ^ a) — - / d 4 a75(l " E «) + 

U c a j \ & a d a 



167T 2 

term of subtraction 



(5.4) 
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where 



(5.5) 



n al = -2[(1 - - a 2 ^ / ]7^[(l - a 2 )i> f - a^] + a^rf + a 2 ^^ f 
+7> x (a l a 3 p* + a 2 a 3 p 2 f + a x oi 2 k 2 ), 

n a2 = a 1 a 2 'fplpf + «?p-^/7 M + «2P/7^/ + ^ 1 a 2 p l i i> f ^p' i p l f , ( 5 -6) 

"l + 20 + 75 \ 



n a3 = [4(1 - ai )pf - 4a 2 p>} - 2mY + + arfj-f] ( + 2 + 7 ' ) 
+ [4(1 - a 2 )pj - 4a lP ? + + aj f y 1 ' + ^ 



(5.7) 



n fl 4 = («iA + (x 2 i> s )i> s Y{^\i>i + ot 2 i> 



1 + 2^-75 
2 



+ [m(«i^ + a 2 ^/)7 M («i + «2^/) + {oL\i>i + a 2 ^/)7 M A(«iA + «2?$/)] 



l + 2fl + 7 5 



n a3 = [4a 3 PZ - 2mY + (a t + a 2 )f o7 M ] ( 1 + ^ + 75 ) 
+ [4«3P<? + 7>i + a 2 )# , o/ 1 



™a4 = 

+QVfo 



( ai + a 2 ) 2 {(fi 2 P o7 M + 2m? ?o - mfiV) ( 3 ) 
1 • 20 - - .7 



n a5 = 2(1 + 20) Jtf - ( 1 + ^ + 7s ) . (5.18) 



(5.8) 

^/l + 2# + 75\ / 1 + 20-75 \ 

n a5 = (2pf - m 7 ^) J+2p^ J, (5.9) 

c a = aia 3 p 2 + a 2 o;3P/ + «i« 2 A; 2 - («i + a 2 )m 2 - a 3 M 2 , (5.10) 

d a = ai(a 3 + «4)p 2 + ct 2 (a 3 + «4)p 2 + aia 2 A; 2 — (an + a 2 )m 2 — a 3 M 2 . (5.11) 

To obtain the term of subtraction in (B.4), we set pi = Pj = Pq with Pq = fl 2 in the 
integrals of (B.4). We get, at the point of substraction, 

c a = («i + a 2 )tt3^ 2 — («i + ct 2 )m 2 — ai 3 M 2 , (5.12) 

d a = (oil + a: 2 )(a 3 + a 4 )fi 2 — (aii + a 2 )m 2 — a 3 M 2 , (5.13) 

n al = -Aa 2 3 P ^ + (1 + a 2 )^V, (5.14) 



n a2 = (ai + a 2 ) 2 ^V, (B.15) 

(5.16) 



(5.17) 
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and 

f oYfo = 2P T o - fiV- (5-19) 
The term of subtraction in (B.4) is proportional to 7 M ^ ^ ^ 5 ^ with the constant of 
proportionality equal to such a value that a + defined in (5.26) is equal to zero at Pi — Pf — Pq. 
Thus, we find, using (B.12)-(B.18), that the term of subtraction in (B.4) is, 

(1 + 19 + 29 2 ) 

g 2 r{l + 2 £ + 75 ) r /■ 3 2 (1 + a ' )fi2 " + + ai + a2)m2 

— {y rf " 2» ( ai + a2 ) a3 W - ( ai + a 2 )m 2 - a 3 M 2 

(1 + 26 + 26 2 ) ^ + a2)(Q2 _ ^ _ a3M 2] Q 2 _ g(1 + g )m 2 Q 2 
-y rf 4 a5(l-^a) p- - - ^ )fi2 _ - ^ 2 _ ^ M2]2 | 

(5.20) 

Next we derive the asymptotic form of Ty^]p in the limit p — > oo with A; fixed. We shall 
only keep the part of this asymptotic form which is relevant to the anomaly amplitude. The 
terms that are proportional to m, i.e. lines 2 and 3 in (B.4) do not therefore, contribute 
to the anomaly amplitude. Let us first observe that in J 2 as given by (4.5), Ty^(pf,pi) 
is always inserted between two propagators. Take, for example, the first trace in (4.5b). 
Let us insert for T p in this trace a next-lowest order term. Then we need to keep only the 
lowest-order terms for S and T p . We shall drop m in S, as the limit p — > oo will be eventually 
taken and one has only to worry about the linearly divergent piece. Thus the first trace in 
(4.5b) is approximately 
i 

We shall keep only the terms in the trace which are linearly proportional to 75, as other 
terms cancel by an extension of Furry's theorem. Since the anomaly amplitude is finite, we 
anticipate the trace in (B.21) be asymptotically proportional to p. (Individual terms may 
actually have factors of hip 2 , but such factors cancel one another.) By counting dimensions, 
we see that terms in the trace in (B.21) must be, asymptotically, proportional to one factor 
of k, where k is an external momentum, in addition to a factor of p. Let us now make 
approximations for the numerators n a \ and n a2 by first moving f { to the right and j>j to the 
left. For example, for n a \ in (B.5), a term i>fi^j)j will be expressed as 

2tff f - 2{ Vl ■ p f )Y + 2p% - tffk. (5.22) 
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r ip + k^p-ki) in p (p + fa) 



(B.21) 



Let us imagine inserting (B.22) into (B.21), with ji,p,i and pf identified with u,p — k\ and 
p + k s , respectively. The third term in (B.22) gives a term proportion to: 



p)Tr 757 W + faW ~ fa) = -P^r 



which is proportional to two powers of k. This means that a power of p has been lost (we 
remind the reader that there is a factor 1/p 2 coming from the denominator c a .) We shall 
therefore drop the third term in (B.22), which does not contain a factor 7^. We also drop 
the first term in (B.22) for the same reason. Thus we get 



Similarly, 



and 



Consequently 



Similarly, 



(^-23) 

i> f ^i> f « -pV (5.24) 

icfii « ~pV (5-25) 
^ / 7 M ^2pV+^/7 M A, (5.26) 
7^/~2pV+^/7^- (5.27) 

n ol » 7 M (2 + « 3 + «3)P 2 +^/7"A(l + «3)- (5.28) 



n a2 « 2 7 ^( P 2 ) 2 (q ;i + a 2 ) 2 +^ /7 ^ 2 («i + « 2 ) 2 - (5.29) 

One must resist the temptation of dropping the last three terms in (B.10) and those in 
(B.ll), as the resulting integrals in (B.4) become divergent. This divergence means that the 
first integral in (B.4) has a lnp 2 term, as the divergence at a 3 = is logarithmic. Similar 
considerations apply for the second integral in (B.4). It is legitimate, however, to drop the 
term a 3 M 2 in (B.10) and (B.ll), and replace pf and p 2 by p 2 . We obtain 

c a ~ ( a l + «2)tt3P 2 + ttitt2^ 2 — («i + ct2)m 2 (5.30) 

and 

da ~ («1 + «2)(«3 + «4)P 2 + «iQ;2^ 2 ~ («1 + a 2 )m 2 . (-B.31) 
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This approximation is allowed because a 3 M 2 is by an order of p 2 smaller than (ai + a2)a 3 p 2 , 
as (cti + CK2) is the order of unity in the region of integration which contributes to the 
asymptotic form. 

Next we make the following change of variables for the second integral in (B.4): 



a 3 
a 3 

«4 



«4 



P, 

px, 

p(l-x). 



(5.32) 



It follows from (B.32) that 



da 3 da^ = pdpdx. 

With the approximation (B.31), the integrands of the second integral in (B.4) is independent 
of x. Thus the intgration over x simply yields unity. We shall rename p as a$, and obtain 
(B.4) as 

.1 . p2 7 M(2_ a3 + a 2) +3 j /7 ^. 



p/x(o)_ 
1 Vipip 



(Pf,Pl) 



9 
16tt 2 



d 3 aS(l - £ 



1 — o^c^p 2 + aia 2 k 2 
(2p 2 Y + ^/7 M A)p 2 a 3 (l - a 3 ) 2 
[(1 — a 3 )o;3p 2 + aia 2 /c 2 — (cki + a 2 )m 2 ] 2 
—term of subtraction. 



(ai + a^m 2 
l + 2^ + 7 5 



(5.33) 



Eq. (B.33) can be rewritten as 



167T 2 
\a1a2k 2 



a) 



p 2 ^a 3 (l + a 3 ) + a 3 $ f ^A 



.(1 — o^c^p 2 + aia^fc 2 — («i + 0.2)™? 
(a 1 + a 2 )m 2 ](l-a 3 )(2p 2 Y+p'fl fl A)] (1 + 29 + 7s 



[(1 — c^o^p 2 + ai« 2 ^ 2 — (aii + ct 2 )'m 2 } 2 
—term of subtraction. 

(5.34) 

We can now drop aia 2 k 2 _ + oi^rn 2 in the denominator of the first term of (B.34), the 
resulting integral being convergent. This cannot be done, however, for the second term of 
(B.34), as the resulting integral remains linearly divergent at a 3 = 0. Instead, this term is 
asymptotically given by the contribution at a 3 = 0, and we may approximate 1 — a 3 by 1 
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and carry out the integration. We get 



2 



J da5(l 



— ^ 5 ^) — term of subtraction (5.35) 



2 



= ^(^7^ + tyfYA/p 2 ) ( 1 + ^ + 75 ) 3 " term of subtraction, 

which is the part of the asymptotic form of contributing to the anolamy amplitude. 

The term of subtraction in (B.35) is the coefficient of the term of subtraction (B.20) with 
a 3 replaced by (1 — a), which depends on the masses m and M. The first term in (B.35) 
comes from the first term in (B.l) which, as we have mentioned, is a convergent integral, 
and has no mass dependences. There is no \np 2 term, as terms of this order cancel. 

(2) Diagram (6b) 

Diagram (6b) gives the following renormalized amplitude to : 

~ lf I C + 2 2 + 75 ) fo + tf-rn^ ~ t6rm ° f Subtraction - 

(5.36a) 

One has to be careful about the term of subtraction here. If we want to make a subtraction 

1 — 7 5 

according to the choice a + = 1, then the divergent part that is proportional to 7 M — - — is 
only partially cancelled by the term of subtraction introduced here. The remaining diver- 
gence will be cancelled by the divergence of the same form in amplitudes corresponding to 
diagrams (6f) and (6g). The term of subtraction, as we introduce here, takes the following 
form 

ip9 ^ N f d A q q 2 + 2a 2 n 2 /l + 20 + 75 



J d 3 ao"(l J 



(1 + 6) J v ^ 'J (2tt) 4 [q 2 + a 2 )a 3 n 2 - (a x + a 2 )m 2 f 

J d 3 a5(l - ]T 



1" 



5.366) 



a )i — ; — — ? — ; — r"?7" 



167T 2 J ^ («! + a 2 )a^VL 2 — («! + a 2 )m 2 



In the asymptotic limit, i.e. p 2 — > oo, we may drop m in the numerator of the integral 
in (B.36a). One then introduces Feynman parameters as usual and make a shift in q, i.e. 
q — > q — a±pi — a 2 p/. The first term in (B.36a) that is relevant to the anomaly amplitude 
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then reads 



d 4 q 2[(1 - a 2 )j> f - axAhKl - - « 2 ^/] - gV 



(27r) 4 [g 2 + aia^pj + a 2 a 3 p 2 + «ia 2 ^ 2 - («i + a 2 )m 2 ] ; 
l + 20 - 75 ^ 



(5.37) 

From the arguments following (B.22), it is easy to see that only the terms $7^ and ^7^ 
in (B.36) contribute to the anomaly amplitude. Using (B.23)-(B.27), the first term in the 
numerator of (B.37) gives 

pyt_%\(l + 29-^ 



-r 



+ 



(5.38) 



2 p 2 /V 2 

while the second term in the numerator of (B.37), together with the first line in (B.36b) and 
part of the subtraction term in (B.45) as will be discussed below in section (5) gives 



16tt 2 



J d 3 a5{l-Y<®) {in 



«3P 



a 3 il 2 — m 2 



r 



l + 20-7 5 



+ (i±^)ln 



+ 



1 + 
1 



a.\ + a 2 )a 3 Q 2 — (cti + a 2 )m 



2 -1 



(ati + a 2 )a 3 fi 2 — 2ax/i 2 — a 3 m 2 



75 



0o^ 2 - (1 + 9)m 



1 + 0/ (cti + a 2 )a 3 f2 2 — («i + a 2 )m 2 



1 + 29 + 7, 



(5.39) 

The contribution of diagram (6b) to relevant to the anomaly amplitude is therefore 

equal to 



-f 

16/T' 



+ 



2 r 



1 

2 



l + 20-7 5 



1 d 3 ao"(l -Ja) 



In 



a 3 p 



a 3 £7 2 — m 5 



7" 



1 + 20 



75 



+ 



1+0 
1 



2 n 



(ai + ct 2 )a: 3 lr - («i + a 2 )m 
«i + a 2 )a: 3 fi 2 — 2ai/i 2 — a 3 m 2 



7" 



75 



0a 2 fi 2 



;i + 0)m 2 



1 + 0/ (ai + a 2 )a 3 f2 2 — (ai + a 2 )m 2 



1 + 20 



(5.40) 



(3) Diagram (6c) 

Diagram (6c) gives the following renormalized amplitude to Vt 



1 tf + * + m / l + 26> + 75 >| ti + 4 + m subtraction 
(2tt) 4 g 2 - 2 W 2 (p f + q) 2 - m 27 I 2 J + g) 2 - m 2 term °* subtractlon - 

(5.41a) 
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Analogous to the discussion for (6b), the term of subtraction in (B.41a) is found to be 



1 + 6) J V ^ J J (2n 



16tt 2 



| d 3 ao"(l - J] 



(2tt) 4 [g 2 + («i + a 2 )a 3 Q 2 - (e*i + a 2 )m 2 - 2a 3 ^ 2 } 

/ 1 + 20 + 75 
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(cti + a 2 )a! 3 n 2 — (ai + a 2 )m 2 — 2a 3 /i 2 



(5.416) 



In a similar fashion, the contribution of diagram (6c) is found to be 



7 M 1 + 2^-75 



+ 



/ 2 



16tt 2 



y d 3 ao"(l - £»{ln 



(ai + a 2 )a 3 Vt 2 — (cm + d 2 )m 2 — 2a 3 fx 2 



/ 1 + 2^-75 



+ 



+ 



1 + 20^ ^ (ai + a 2 )a: 3 £r - («i + a 2 )m 2 - 2a 3 ji 2 ^^(1 - 75 



1 + 
1 



(«i + a 2 )a: 3 f2 2 — 2a 1 /i 2 — a 3 m 2 
0a 2 ft 2 - (1 + 9)m 2 



1 + 6/ (ati + a 2 )a 3 Vt 2 — (ati + a 2 )m 2 — 2a 3 /z : 



r 



l + 20 + 7 5 



(5.42) 



(4) Diagrams (6d) and (6e) 



Diagram (6d) and diagram (6e) give the following amplitude to valid to the next- 

lowest order: 



2ig 2 M 



A r ti + i + m ,(1 + 20 + 75) ,(1 + 20 + 75) j)f-i + m 
7 ^ ^7 



(2tt 



( Pl + q) 2 -m 2 ' 2 1 2 (p f - qf 



m 



gla , _ qf , q u /g 2 x 

n ttk — -7-, o — term or subtraction. 

q 2 — M 2 (k + q) 2 - 2/i 2 



(5.43a) 



The term of subtraction in (B.43a) is evaluated to be 

-2g 2 Mm f r 3 ^ 2j(l~Ea) 

167r 2 I 7 a^ 2 — Q! 1 (l — a^)!] 2 + 2a 2 fi 2 + a 3 M 2 

<$(!-£) a) 1 u/ 1 + 20 + 75\ ^ ' ' 



+ ( d A a- 
J ( 



aim 2 - ai(l - a^Q 2 + 2a 2 /i 2 + a 3 M 2 J V 2 
Because of the factor of M outside the integral sign, we conclude by counting dimension that 
the first term in (B.43a) does not contribute to the anomaly amplitude. Thus the expression 
in (B.43a) is replaced by the negative of (B.43b). 

(5) Diagrams (6f) and (6g) 



38 



Diagrams (6f ) and (6g) give the following renormalized amplitude to 



A , n _ , ,_\u\ A + 4-m 1 1 



^ pi + 4 + m 1 



(pi + q) 2 — m 2 q 2 — 2\i 2 (k + q) 2 



7 5 — term of subtraction. 



(BAA) 



(Pi + <?) 2 — m2 (f (k + q) 2 — 2/j, 2 . 
where, in the first term of the integrand, there is a matrix 75 which we have moved to the 
rightest position. The integral in (B.44) has a logarithmic ultraviolet divergence, which is 
not entirely cancelled by the term of subtraction in (B.44). This is because the former is 
proportional to 7^75, which can be written as: 

7^(1 + 75) -7^(1 -75). 

The term of subtraction that we introduce here only cancels the divergent part of the form 
7 M ( — - — J. The uncancelled divergence, in the form of — ^ 5 )' * s cancene d by the 

divergence of the same form in amplitudes corresponding to diagrams (6b) and (6c) that still 
remains as we have discussed above. An amplitude of such a form is ultraviolet finite and 
requires no subtraction, as we have concluded in Sec. 2. 

Let us single out the terms in (B.44) which contribute to the anomaly amplitude. First 
of all, we drop m in the numerators in (B.44). Then the relevant part in (B.44) is 



"7 7^< 2 « + ^ 



-term of subtraction. 



1 1 

+ 



q 2 - 2/i 2 (k + q) 2 q 2 (k + q) 2 - 2/j, 2 



75 



(5.45) 



Introducing Feynman parameters and calling 

q =q + 5q, 

5q = -a 3 pi - a 2 k, 

we find that, since an amplitude of the form of ± 75) does not contribute to the 

anomaly amplitude, as we have discussed in the paragraph following (B.22), we may drop 
all terms in (B.45) involving 5q. The subtraction term is given by 

*f 2 f ,3 c n _ v x V / l + 2g + 75 \ 

(1 + 6) J 1 ^ a> J (2tt) 4 [q 2 + + a 2 )a 3 Q 2 - 2a lf i 2 - a^m 2 ]* 1 V 2 ) ' 

(5.46) 
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We break up the last factor in (B.46) into two terms: 

l+2i+Ji = (1 + 0h5 + (1 + 20) 



1 - 75 



(5.47) 



Then the subtraction term in (B.46) becomes a sum of two terms. The first term combines 
with (B.45) to give the following contribution to the anomaly amplitude : 



2p_ 
16tt 2 



a 3 m 2 — («i + a 2 )a 3 fl 2 + 2ai/i 



2n 



7^75- 



(5.48) 



-(«i + a 2 )a 3 p 2 

The remaining term combines with the subtraction terms in (B.36b) of diagram (6b) and in 
(B.41b) of diagram (6c) to give a finite contribution to the anomaly amplitude as given in 
(B.40) and (B.42) above. 
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Appendix C 



In this appendix, we discuss the lowest-order radiative corrections of the fermion propa- 
gator, with special emphasis on their contributions to the anomaly amplitude. The relevant 
diagrams are illustrated in Fig. 11. 

(1) Diagram (11a) 

The renormalized amplitude corresponding to diagram 11a is: 

T (r) (p) = if 2 ( f 1 t-a t~~ — m n n — n — term of subtraction. (CM) 
^ a KV> J J (2tt) 4 (p - q) 2 - m 2 q 2 - 2/i 2 V ; 

Introducing Feynman parameters, we get 

rw \ • r2 f 1 t f d A Q (l — a)4> + m „ . 

£«(p) = «/ J q da J Wr * [? + a (l- a )p>- am *-2{l-aW]* " ^ ° f Subtractlon ' 

(C7.2) 

where 

q = q — ap. 

The integral in (C.2) is logarithmically divergent. After we make a subtraction it is equal 
to: 

Ma(p 2 )+mB a (p 2 ), (C.3) 



where 



and 



-/ 2 _ Nl _f a 2 m 2 + 2(1 - a)^ 



Mp2) = m^) I da{1 - a) 111 



(167T 2 ) Jo 

Baip2) = a&) Jo da 111 



— a(l — a)p 2 + am 2 + 2(1 — a)/i 



f 2 f 1 J , r a 2 m 2 + 2(1 - a)/i 2 



(16tt 2 ) 

The subtraction is so chosen that 



— a(l — a)p 2 + am 2 + 2(1 — a)ji 2 



(C4) 



(C.5) 



A a (m 2 ) = 5 a (m 2 ) = 0. 



Expression (C.3) is not yet the answer for J2a\ the subtraction for which is given by 
(2.14). Following the prescription of (2.14), we get 

= Ma(p 2 ) + mB a (p 2 ) - 2m 2 K(m 2 ) + # a (m 2 )]0 - m). (C.6) 
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From (C.4) and (C.5), we get 



-f 



Jo 



da- 



ail — a)' 



and 



(167T 2 ) Jo a 2 m 2 + 2(1 — a)// 2 ' 
■P f 1 a(l-a) 



/' 

JO 



' a 2 m 2 + 2(1 — a)/x 2 ' 
It is then straightforward to find that, as p — > oo, 



El r) (p) 



if 

16tt 2 



i (4 - 4a + a 2 )m 2 - (2 - a)^ 2 



-ln(— p 2 /m 2 )+ I daa- 

2 y F 1 'Jo a 2 m 2 + 2fi 2 (l -a) 

(2) Diagram (lib) 



(C.7) 



The renormalized amplitude corresponding to diagram (lib) is 



Er ) (p)=</ 2 



<i 4 g ^ 



m 1 



(27r) 4 (p — g) 2 — m? q 2 
Comparing (C.8) with (C.l), we find that 



term of subtraction. 



£(-)(p) = p'Aaip 2 ) - mB a (p 2 ) - 2m 2 [A' a (m 2 ) - B a (m 2 )}(^ - m), 
with jj, in A a and _B a set to zero. Thus we get, for p — > oo, 



E1 p) (p) 



-ln(-p 2 /m 2 ) - - 



(16vr 2 ) 
(3) Diagram (11c) 



(C.8) 



(C.9) 



(CIO) 



Finally, we turn to diagram (11c). The renormalized amplitude in the Landau gauge 
corresponding to this diagram is 

A (3j - 2^ - flrf /g ») 1 + 2g + + (1 + ^ + 3^(1 + fl) 
^ 7 (2tt) 4 [(p-g) 2 -m 2 ](g 2 -M 2 ) 
—term of subtraction. 

(C7.ll) 

Introducing Feynman parameters, making a shift of the variable, and carrying out the inte- 
gration, we get 



A c {p 2 )j> 



1 + 29 + 29 2 + (1 + 2#) 75 



+ mB c (p 2 )6(l + 6) 



(C.12) 



—term of subtraction, 
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where 



c{P 1 " Jo J (2tt) 4 [q 2 + a(l - a)p 2 - am 2 - (1 - a)M 2 } 2 



q — 2a\p 



(C7.13) 



(2tt) 4 [q 2 + ai(l - «i)p 2 - aim 2 - a 2 M 2 ] 3 ' 
While each of the two integrals in (C.13) is logarithmically divergent, the sum of the two 
integrals is ultraviolet finite. Carrying out the integration over q we get, after some algebra, 

1 f 1 



Mp 2 ) 



16tt 2 



da 



(2a — I — a 



2 2 

p — m 

M 2 



) ln[-a(l - a)p 2 + am 2 + (1 - a)M 2 } 



+a 



p 2 — m 2 
M 2 



ln[— a(l — a)p 2 + am 2 ] — a 



(C.14) 



From (C.14), it is straightforward to obtain 

A' c (m 2 ) 



1 r i « 2 ("^ + 2«) 
da- 



167T 2 Jo a 2 m 2 + (1 — a)M 2 
In a similar way, B c (p 2 ) and B' c (m 2 ) are given by 

r d A q 3 



(C.15) 



(2tt) 4 [g 2 + a(l - a)p 2 - am 2 - (1 - a)M 2 ] : 



— term of subtraction 



167T 2 JO 



da In 



l^(m 2 ) 



-a(l - a)p 2 + am 2 + (1 - a)M 2 
a 2 m 2 + (1 - a)M 2 . 

3 r 1 , a(l -a) 



f 

Jo 



da- 



167r 2 Jo a 2 m 2 + (1 - a)M 2 ' 
From (C.14), (C.15), (C.16) and (C.17), we get, following the prescription of (2.14), 



(C.16) 
(C17) 



E1 p) (p) = 



, , -ail - a)p 2 + am 2 + (1 - a)M 2 

cto(2a - 1) In t—, 

' a 2 m 2 + (1 - a)M 2 

-a(\ - a)p 2 + am 2 1 1 + 20 + 26 2 + (1 + 2#) 75 



16tt 2 

Qo 2 - m 2 ) 

" M 2 n -«(l-«)p 2 + «m 2 + (l-a)M 2 . 

^2 m 2 i a 2 (-^ + 2a)(l + 2# + 2# 2 ) + 6a(l-a)#(l + #) 

+ l6Wo ^ a 2 m 2 + (1 _ a)M 2 ^~ m ) 

3g 2 m n ^ m f 1 , , [ a 2 m 2 + (l-a)M 2 



16tt 2 



0(1 + 0) / da In 

JO 



— a(l — a)p 2 + am 2 + (1 — a)M 2 



(C.18) 

The part of relevant for the anomaly amplitude is obtained from (C.18) by taking 
the limit p — > oo. We find, after some algebra, that the part in which contributes to 
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the anomaly amplitude is 

C da—- i — - ( [2a 2 (2 - a)m 2 + (1 - a 2 )M 2 ] (1 + 29 + 29 2 ) 

32n 2 Jo a 2 m 2 + (l-a)M 2 V V ; V ; JV ; (^.19) 

+ [a 2 (2a - l)m 2 + (1 - a 2 )M 2 ](l + 26)^ 
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Appendix D 



In this appendix we evaluate the ratio Z^(m) / Z V ^(Q, , , 0). This quantity will be 
evaluated by the use of (2.27). 

The values of c r (Q 2 ), d r (Q 2 ) and a r (Q 2 ) have been given in Appendices B and C, therefore 
we only need to calculate G + (Q 2 , Q 2 , 0). After some calculations, we find that 



(«i + a 2 )m 2 - (cti + a 2 )a 3 fi 2 + 2a 3 /x 



2 



1 2/i a 2 — a\ 



Oi\ + a 2 )m 2 — («i + a 2 ) 2 f2 2 ot.^m 2 + 2a 2 /i 2 — («i + a 2 )a: 3 f2 2 

mg 2 /■ 3 ^ f 2 

16^7 1 ~ ^ aj I c^m 2 - aAl - a^fl 2 + 2a 2 pi 2 + a 3 M 2 

1 1 

2 («! + a 2 )m 2 + a 3 M 2 — (ai + a 2 )a: 3 n 2 

W [ ,4 r/-, \ - xf 2 

16tt 2 7 1 Z - aj U 1 m 2 -a 1 (l-a 1 )fi 2 + 2a 2 /i 2 + a 3 M 2 
+ 



2 («i + a 2 )m 2 + a 3 M 2 — (a± + a 2 )a 3 f2 2 



(D.l) 
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Appendix E 



In this section, we shall calculate the anomaly amplitude for Ji and J 2 . The calculation 
of the volume integral J\ will be presented in Section (a), while the calculation of the surface 
integral J2 will be presented in Section (b). 
(a) The Volume Integral J\ 

The Feynman integral corresponding to each of the diagrams in Fig. 7 is quadratically 
divergent by power counting. However, the trace in the integrand of each integral vanishes 
if we set all of the momenta of the external vector mesons to zero. This means that each 
term in the trace contains at least one factor of the external momenta and hence the trace 
blows up no faster than the cubic power of the internal momenta. Consequently, each of the 
Feynman integrals in question is linearly divergent. As we have discussed in Section 3B as 
well as in Section 5, the sum of these integrals are convergent by the requirement of Bose 
symmetry. 

We shall first carry out the integration over q. We shall do this for the sum of anomaly 
amplitudes corresponding to diagrams (7a) and (7b). As we can conclude from power count- 
ing, this subintegration over q is convergent. Next we carry out the integration over p. As 
we have discussed in detail in Section 5, this integration is also convergent as we sum over 
the contributions from all the diagrams. The anomaly amplitudes corresponding to dia- 
grams (7c) and (7d) are equal to those corresponding to diagrams (7a) and (7b) with the 
interchange of 2 <-> 3. We shall show that this sum is equal to 

J i = -j^y^^'hahAii + ef - 9% (e.i) 

a constant which depends on the renormalized coupling constant g, and not on the masses 
of the particles in the theory. This is because the contributions to these anomaly amplitudes 
come from a region in which the internal momenta are infinitely large, hence the particle 
masses which appear in the integrand can be dropped. We shall next present the calculation 
of these amplitudes. Readers who are not interested in this calculation are advised to skip 
the rest of the subsection. 

Referring to diagrams (7a) and (7b), we find the sum of amplitudes from these two 
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diagrams as 

where 

N A = Tr 



f d A p d A q 






1 (27r) 4 (27r) 4 







(E.2) 



„ l + 2fl + 75 ^ 1 + 20 + 75 , v ^ pI+20 + 75,, « 

7 o {p + hn o {p + h-in p o O-fi 



7 U 



CT ,/l + 2^ + 75\ 



J J* 



2 

w 5— 

v Q 



(E.3) 



V A = (p 2 — rn 2 ^j (p + A; 2 ) 2 — m 2 (p + k 2 — q) 2 — m 2 (p — k± — q) 2 — m 2 ( 



(£.4) 



and similarly for N B and Dg. We have dropped the fermion mass m in the numerators of 
the fermion propagators. This is because the terms that involve m in the numerator are 
convergent integrals. Thus, shifting the loop-momenta for the terms which have m in the 
numerator is allowed and hence such terms in M A cancel the corresponding terms in Mb- 
The propagator of the vector meson is in the Landau gauge, which accounts for the last factor 
in (E.3). As noted, the contribution of the left-handed fermion and the right-handed fermion 
can be calculated separately, as all five vertices on the fermion loop are Vi/j^ vertices. Thus 



with 



and 



N 



(i) 



Tr 



i^ v (P' + h)(i>-h-i)i p (p' + h-i)p' 



N 



(2) 



(E.5) 



(E.6a) 



(E.6b) 



where we have averaged over the two directions of the fermion loop. We shall put 

M A = M A l) + M A 2 \ 



where Xl ( } ] { M { A ] 



is obtained from M A by replacing N A with with N^^N A '/q' 

r(2) 



r(2) 



We 



also define M B and M B similarly with reference to diagram (7b). By introducing Feynman 
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parameters and carrying out the integration over q, we get, after some algebra, 

d 4 p 1 



rN 



(i) 



(2tt) 4 A 3 

a 3 a b Tr(^ u l p hh) 



(p - A A ) 2 + rD A A~ 



2r 



(p - A A ) 2 + r^A^ 2 " 
a4a 5 Tr( 7 57"7^2^3) 



(p- A B ) 2 + rD B A 



-2 



2r 



(p — Ab) 2 + r Db/\ 



> . 



(£.7) 



In (E.7), 



D A = Q;iQ;3Q;5/c 2 ai + [a 5 (a 2 + CK4) + «2(«3 + 0:4)] &2 + a 3 [(<*1 + «2)«4 + («2 + «4)«5] &| 

+ [(ai + a 2 + a 3 + a 4 )m 2 + a 5 M 2 ] A, 
A = (ax + a 2 )(a 3 + a 4 ) + a 5 (ax + a 2 + a 3 + a 4 ), 
r = a 3 + a 4 + a 5 ; 



also, 



with 



= Tr 



p A =p- A A , 



and 



Api = ai(a 3 + a 4 + a 5 )k 2 — a 3 a 5 k 3 , 

Ap 2 = a 1 a 5 k 2 + [(a 1 + a 2 )(a 4 + a 5 ) + a 4 a 5 ]k 3 , 

Ap 3 = a^ 5 k 2 — a 3 (a.i + a 2 + « 5 )/c 3 , 

AAy! = -[(« 2 + «5)("3 + "4) + a2«5]^2 - Oi 3 a 5 k 3 . 

In addition, D B is equal to .D^ with cti <-> a 2) ^3 <-> a 4 : 



Cti <-> CK2, CK3 <-> 0:4 



(£8) 



(£9) 



(£.10) 



ivL 1} = Tr 



r r 



(£.11) 
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with 

p B = p - A B , 

Ap[ = ^40:5^3 + + a 5 )(o;3 + a 4 ) + aia 5 ]/c2, 

Ap' 3 = -a 2 tt5^2 + [«3(«i + "2) + a 5 (o;i + a 2 + a 3 )]k s , 

AA B = «4tt5^3 — «2(«3 + «4 + «5)&2- 

(£.12) 

Next we carry out the integration over p in (E.7). Consider the first term of the integrand 
in (E.7). We shall change the variable of integration from p to pa which is defined in (E.10). 
This is to say that, instead of performing a symmetric integration over p, we shall perform 
a symmetric integration over pa- Since N A ^ given by (E.9) is cubic in pa, the corresponding 
integral is linearly divergent, and a term proportion to A^ is generated by the change of 
variables. In addition, we note that the term in N A ^ obtained by setting pa to zero is exactly 
canceled by its counterpart in M$. This is because diagram (7a) is topologically identical 
to diagram (7d), and the amplitudes corresponding to these two diagrams are different only 
because the loop-momenta are chosen differently. But p\ in N A ^ is equal to its counterpart in 
N£\ To see this, let us think of diagram (7a) as an electrical circuit, with a, the resistance 
of line i. Then p 1 for example is the current in line 2. Since diagram (7d) represents the 
same electrical circuit, it has the same currents in its lines. This is why N A ^ at p A = is 
exactly equal to N$ at p D = 0. Thus we shall ignore such a term in N A \ Similarly, we 
shall ignore the term in N B ^ which is independent of pb- After dropping such terms and 
calling both pa and p B as p, the bracket in (E.7) becomes, after some algebra and changes 
of variables, 

p 2 a 5 [2a 3 (a 2 + 05) - 2a 1 (a 4 + a 5 ) - a^r + Q^lA"] c^c^ . „ . 

N{f + ZV/A 2 ) 3 r{p 2 + D a r/A 2 y 1 ' j 

Substituting (E.13) into (E.7) and carrying out the integration over p, and adding to it 
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— M$ which is obtained from — by the interchange of 2 <-> 3, we get 



M 



(i) 



m£ ] + M, 



(i) 



c 



M 



(i) 



d 5 o<5(l-I» 
A 1 " 



0:10:5(0:3 + 04 + 05) — 0305(03 + 04 + 05) + 305(03 — Oi)/n ^ 
05(03 + 04)(o3 + 04 + 205) 1 Ao 5 (o 3 + o 4 



Oi + o 2 + o 5 
03 + o 4 + o 5 



o 3 + o 4 + o 5 



2 03 + o 4 + o 5 



(£.14) 

with the last two terms in the integrand derived from the shifting of the integration variables. 
The integrations over the Feynman parameters are routine and we find after some algebra 
that the right-side of (E.14) is equal to the expression in (E.l). 



Next we calculate Nlf ] - . We get 



Mf - M<?> 



2 y 



d A p d A q 1 
(27r) 4 (27r)V 



N 



(2) 



N 



(2) 



B 



V, 



with Ny given by (E.6b), and with N%' given by a similar expression with reference to 
diagram (7b). 

The ensuing calculations are somewhat simplified if we replace the first ^ in (E.6b) by 



Then we get 



jVi 2) = (p + k 2 ) 2 Tr 



(p + k 2 — q) 2 Tr 



Now 



(p + /c 2 ) 2 = [(p + A; 2 ) 2 - m 2 ] + m 2 . 



(E.15) 



(E.16) 



The first term on the right-side of (E.16) will be used to cancel the same factor in P4. The 
second term in (E.16), m 2 , will be dropped. This term is by two powers of loop-momenta 
less divergent than the first term, and its contributions are cancelled by its counterpart of 
diagram (7d). Next we introduce Feynman parameters, integrating first over q and then 
symmetrically over p. In the latter integration, a change of variable from p to p is made. 
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The part of — + M^) — M$ which comes from symmetric integration over p is 
calculated to be 

2ig r daida 2 da 3 da55(l — J2 a)a\ 

-Tr (757 Yfcfa) 



(167T 2 ) 2 J [(«! + o 2 )o 3 + a 5 (ai + o 2 + o 3 )] 

/« 1 + «2 + « 5 \ , a 1 (a 1 + a 2 )a 5 o 2 o 5 

-3a 3 (Q!i + a 2 ) m ■ H ■ ■ ■ o 2 o 3 . (E.17) 

\ 0:3 + 0:5 / Oi + o 2 + o 5 o 3 + o 5 

It is straightforward to carry out the integration over the Feynman parameters. We find 
that the expression in (E.17) is equal to 

-Ug 2 (lQn 2 )- 2 Tr( l5 YYhh)- (E.18) 

In addition, there is a term which comes from the change of integration variables from p to 
p. This term is equal to the negative of (E.18). Thus 

il4 2) - il4 2) + il4 2) - M§ ] = 0. (£.19) 

Consequently, the sum of anomaly amplitudes corresponding to the diagrams in Fig. 7 is 

[(1 + ef - e 5 ][M A -M B + M C - m d ] = -^L^ g 2 e u P a(T 'k 2a k 3a/ [(i + ef - e 5 ]. (E.20) 

(b) The Surface Term J 2 

In this subsection we evaluate the surface integral J 2 given by (4.5) in the two lowest 
orders of the renormalized coupling constants. 

First we express the unrenormalized functions T v ^ and S in (4.5b) by the corresponding 
renormalized ones. Then N-j is given by (4.5b) with T and S replaced by and S<?\ 
respectively, with an additional overall factor 



Z4m)/ Z v ^(Q 2 ,Q 2 ,0) 



(£•21) 



All of the fermion momenta at the surface of integration of (4.5a) are infinitely large. There- 
fore, we only need to consider the asymptotic forms of T and S in this limit. Referring to 
(2.2) and (2.12), we have, for p 2 — > oo, 

L[S (r) (p)f 1 « -il4{\ -c r + d r ), (E.22a) 

[S (r) (p)f *L w -c r - d r )L, (E.22b) 
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and 

P[S (r) (p)] _1 « -iRfi(l -cr- dr), (£.23a) 

[S ,(r) (p)] _1 « -^(1 - c r + d r )#. (£.236) 

where L and R are defined in (3.16). As to r^^, it is a linear superposition of a number of 
invariant amplitudes each multiplied by one of the following terms: 

u 1± l5 m „1±75 v 1/ 1 i 75 // z,/ 1 ±75 

A „l±75 ftA „l±75 - AI/ 1±75 jl/A'U 1± 7 5 

mA , f A , AA — — , A ^ , 

mP — ^— , PP — ^— , AP — ^— , P • 

(£.24) 

All terms in (E.24) are of the dimension of even powers of momentum. In this way, the 
invariant amplitudes associated with these terms are also of the dimension of even powers 
of momentum, and are functions of p 2 ,p' 2 and A 2 . 

The function T v ^ is dimensionless. Therefore, in the limit of p 2 and p' 2 going to infinity 
with A 2 fixed, the surviving terms in (E.24) are 

7^^4p , , PP'^4p. (£.25) 

1 ± 75 

For example, the coefficient of by dimensional considerations of the order of 

_ 1 ± 7 5 

(P 2 ) -1 as P 2 — > 00 . Since (P 2 )~ 1 mf / y u vanishes as P 2 — > 00, this amplitude can 

2 

be neglected. Thus, all the terms in (E.24) with a factor m can be neglected. Similarly, 

all the terms in (E.24) with a factor A can be neglected. We may also show that the term 
1+75 

fP v — - — in (E.25) does not contribute to J 2 , as it is possible to verify that N 7 vanishes 
asymptotically if we substitute such a term for T v ^ in (4.5b). Thus, the terms in 

(r) 
Vtptp 



contributing to J 2 are the first two terms of (E.25). We shall therefore put in the form 



II + a+ )7"(l + f) 1 -^ + (1 + a-W-^*- + rjl±(l + «)i±* + 

(£.26) 
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Substituting (E.23) and (E.26) into (4.5b), we find after some algebra that, in the two lowest 
orders in the coupling constant 



J 2 = e 



i/paa 



i klcrk2a' 



Z^(m) 



1 2 



'1 + 9) 3 [(1 + 2a + - T -±- + 2c r + 2d r ] 

1 (E.27) 



-9 3 [1 + 2a_ - T — + 2c r - 2d r ] J. 



The asymptotoic forms for a±, c r and d r have been evaluated explicitly in Appendix B and 
C, while the value of Z 1 p(m)/Z v ^(Q 2 ,Q 2 ,0) was evaluated explicitly in Appendix D. The 
calculations in these appendices are straightforward but tedious. Such explicit calculations 
can be avoided if we make use of the Ward-Takahashi identity (2.19) in the limit of p 2 and 
p' 2 going to infinity with A 2 fixed. Let 



= 1 + 2. (E.28) 



Substituting (E.26) into the Ward-Takahashi identity (2.19) and making use of (E.22) and 
(E.23), we can equate the coefficients of A ^ ^~^ 5 and obtain, in the limit of p 2 and p' 2 going 
to infinity with A 2 fixed, 

z + a± - t± + c r ± d r = 0. (£.29) 



Substituting (E.29) into (E.27), we get 

3 3 _ 1 

(£.30) 



(1 + 0)3 _^ + 3( 1 + 0)3 T+ _303 T _ 



In Appendix B, we showed that the only diagram which yields an amplitude r + are diagrams 
(6a), (6b) and (6c). The value of r± are given by 
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Figure Captions 

Figure 1. Lowest-order triangular graphs. The thick straight lines represent fermions, and 
the wavy lines represent vector mesons. Only two of the six diagrams are 
topologically independent. 

Figure 2. A double line represents a 2 -meson. See the caption of Figure 1 for the meanings 
of other lines. 

Figure 3. A vertex with a blob represents a dressed vertex, with all of the radiative corrections 
included. A line with a blob represents a dressed propagator. 

Figure 4. A straight thin line represents a H-meson. For meanings of other lines, see caption 
of Figure 1. 

Figure 5. A diagram in this figure represents an anomaly amplitude. A dotted line in 
these diagrams represents the vector meson of momentum ki, the vertex for 
which is absent as the Ward-Takahashi relation (2.17) has been applied. 

Figure 6. The diagrams for the lowest-order radiative corrections of the Vi/jiJj vertex. 

Figure 7. The anomaly diagrams obtained from those in Figure 5 with a dressed vertex 
replaced by a vertex graph in Figure 6a. 

Figure 8. The anomaly diagrams obtained from those in Figure 5 with a dressed vertex 
replaced by a vertex graph of Figure 6b. 

Figure 9. The set of nonplanar diagrams which contributes to the next order 
triangle anomaly. 

Figure 10. Diagrams depicting the reduction of amplitudes by the use of (4.8). A dot on 
an external vector-meson line represents the multiplication of ik^ to the Vipip 
amplitude in (4.8). 

Figure 11. Same as Figure 10 but with mesons 2 and 3 interchanged. 

Figure 12. The anomaly diagrams obtained from those in Figure 5 with a dressed vertex 
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replaced by a vertex graph in Figures (6f) and (6g). 

Figure 13. The anomaly diagrams obtained from those in Figure 9 which contribute to 
the next order triangle anomaly amplitudes. 

Figure 14. Diagrams representing three anomaly amplitudes which cancel one another. 

Figure 15. Lowest-order radiative corrections for the fermion propagator. 
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